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During  this  contract  period,  substantial  progress  was  made 
in  the  development  of  a  comprehensive  model  for  simulating  hot 
fragment  conductive  ignition  (HFCI)  of  low- vulnerability  ammuni¬ 
tion  (LOVA)  propellants.  A  numerical  code  was  also  developed  to 
solve  the  HFCI  model.  The  work  has  resulted  in  three  major  pub¬ 
lications,  included  herein  as  Attachments  1,  2  and  3. 
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(acceptance  for  publication  included  as  Attachment  4) 
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Attachment  3:  Hsieh,  K.  C.,  Hsieh,  W.  H. ,  Kuo,  K.  K. ,  and 
Miller,  M.  S.,  "Validation  of  a  Theoretical  Model  for  Hot  Frag¬ 
ment  Conductive  Ignition  Processes  of  LOVA  Propellants,"  Proceed¬ 
ings  from  the  10th  International  Symposium  on  Ballistics  ( San 
Diego,  California,  October  27  -  29,  1987),  Vol .  1,  Section  2. 

The  model  and  code  developed  under  this  contract  can  be 
applied  in  studying  the  vulnerability  of  various  propellants 
heated  by  spall  fragments. 
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ABSTRACT 

A  comprehensive  theoretical  model  has  been 
formulated  for  studying  the  degree  of  vulnerability  of 
various  solid  propellants  being  heated  by  hot  spall 
fragments.  The  model  simulates  the  hot  fragment 
conductive  ignition  (HFCI)  processes  caused  by  direct 
contact  of  hot  inert  particles  with  solid  propellant 
samples.  The  model  describes  the  heat  transfer  and 
displacement  of  the  hot  particle,  the  generation  of 
the  melt  (or  foam)  layer  caused  by  the  liquefaction, 
pyrolysis,  and  decomposition  of  the  propellant,  and 
the  regression  of  the  propellant  as  well  as  the  time 
variation  of  its  temperature  distributions. 

To  partially  validate  the  theoretical  model  in 
the  absence  of  the  necessary  chemical  kinetic  data,  an 
ice  melting  and  evaporation  experiment  was  designed 
and  conducted.  These  experiments  provide  features  of 
the  conductive  heating,  melting,  and  evaporating 
processes.  Calculated  results  compare  well  with 
experimental  data  in  temperature-time  traces,  spall 
particle  sinking  velocity,  and  displacement. 

NOMENCLATURE 

A  Arrhenius  frequency  factor  of  chemical 

reaction  of  melt  or  foam  layer  in  Region  1, 
kg/m^-s 

A-j  Interface  area  between  liquid  melt  and  gas 

bubbles,  m2 

Cp  constant-pressure  specific  heat,  J/kg-K 

Cpr  specific  heat  of  LOVA  propellant,  J/kg-K 

Cs  specific  heat  of  spall  particle,  J/kg-K 

Ob  averaged  diameter  of  bubbles  generated  from 

gasification  process,  m 

Eg  activation  energy  of  chemical  reaction  of 

melt  or  foam  layer,  J/mole 
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convective  heat  transfer  coefficient  on  top 
and  lateral  surfaces  of  spall  particle, 
W/tn2-K 

convective  heat  transfer  coefficient  on  top 
surface  of  LOVA  propellant,  W/m^-K 
heat  of  formation,  J/kg 
latent  heat  of  liquid  melt,  J/kg 
instantaneous  height  of  melt  or  foam  layer 
in  Region  2,  m 

instantaneous  distance  traveled  by  spall 
particle,  m 

height  of  melt  or  foam  layer  in  Region  1,  m 
height  of  the  spall  particle,  m 
rate  of  conversion  of  reactant  (R)  into 
product  (P),  kg/s 

rate  of  mass  generation  of  reactant  species 
from  the  base  of  spall  particle,  kg/s 
net  rate  of  production  of  gaseous  mass  in 
Region  1  of  foam  layer,  kg/s 
order  of  chemical  reaction  in  gas  bubbles 
pressure,  N/m^ 

critical  heat  flux  for  producing  gaseous 

bubbles  from  melt  of  LOVA  propellant  on 

surface  of  spall  particle,  W/m^ 

heat  flux  on  submerged  lateral  surface  of 

spall  particle,  W/m^ 

heat  flux  from  foam  layer  to  LOVA 

propellant,  W/m^ 

net  heat  generation  in  LOVA  propellant, 

W/m^ 

heat  flux  transferred  to  LOVA  propellant 
from  melt  or  foam  layer  in  Region  1,  W/m^ 
net  radiant  heat  flux,  W/m^ 
heat  flux  at  base  of  spall  particle,  W/m^ 
radial  coordinate,  m 

Universal  gas  constant,  8314.4  J/kg-mole-K 
outer  radius  of  melt  or  foam  layer  in 
Region  2,  m 

radius  of  spall  particle,  m 
boiling  temperature  of  liquid  melt,  K 
gas  temperature,  K 

bulk  temperature  of  melt  or  foam  layer,  K 
melting  temperature  of  LOVA  propellant,  K 
temperature  of  LOVA  propellant,  K 
average  surface  temperature  of  LOVA 
propellant  exposed  to  air,  K 
temperature  of  spall  particle,  K 
average  axial  velocity  of  bubbles  in  Region 
1  of  foam  layer,  m/s 

radial  velocity  of  melt  or  foam  layer  in 
Region  1,  m/s 

axial  velocity  of  boundary  surface  of  melt 
or  foam  layer  in  Region  1,  m/s 


''mr 


outward  radial  velocity  of  foam  layer  in 
Region  2,  m/s 

sinking  velocity  of  spall  particle,  m/s 
axial  distance  above  base  of  spall 
particle,  m 

axial  distance  below  instantaneous  surface 
of  propellant  under  foam  layer,  m 


Greek  Symbols 


a 

B 

A 

P 

es 

0 

Os 
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thermal  diffusivity,  m^/s 
bubble  contact  angle,  degree 
thermal  conductivity,  W/m-K 
density,  kg/nr 

average  porosity  (void  fraction)  of  foam 
layer  in  Region  1 
emissivity  of  spall  particle 
Stefan-Boltzmann  constant,  ■  5.6696  x  10“® 
W/m2-K^ 

surface  tension  of  liquid  melt  of  LOVA 

propellant,  N/m 

dynamic  viscosity,  kg/m-s 

average  void  fraction  of  spall  base  surface 

in  contact  with  gas  bubbles 


Subscripts 

g 

i 

1 

melt 

pr 
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gas-phase 

initial  condition 

liquid-phase 

liquid  melt 

propellant 

room  condition 

region  1  of  foam  layer 
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*  T*  -  T  -  Tref 

INTRODUCTION 


Ignition  of  solid  gun  propellants  in  go/no-go 
tests  which  employ  a  hot  metallic  element  of  well- 
defined  geometry  as  the  source  of  energy  has  been  a 
conventional  method  of  determining  the  vulnerability 
of  propellants  (Gol 'dshleger,  1973),  especially  in  the 
Low  Vulnerability  Ammunition  (LOVA)  development 
program  (Wise  et  al.,  1980;  Law  and  Rocchio,  1981;  and 
Wise  and  Rocchio,  1981).  Hot  fragment  conductive 
ignition  (HFCI)  studies  help  to  determine  the 
survivability  of  weapon  systems  (such  as  tanks,  ships, 
etc.)  containing  stowed  ammunition  (Gol 'dshleger  et 
al.,  1973;  Wise  et  al.,  1980;  Law  and  Rocchio,  1981; 
and  Wise  and  Rocchio,  1981).  Ignition  of  propellant 


charges  by  hot  metallic  spall  fragments  can  be 
generated  through  such  threats  as  shaped  charges  and 
kinetic  energy  penetrators. 

Review  of  HFCI  processes  indicates  that  both 
experimental  and  theoretical  investigations  have  been 
conducted.  The  experiments  generally  employ  hot  steel 
balls  of  known  diameters  and  temperatures  which  are 
suddenly  brought  into  contact  with  the  propellant 
(Gol 'dshleger  et  al.,  1973a;  Wise  et  al.,  1980;  Law 
and  Rocchio,  1981;  and  Wise  and  Rocchio,  1981). 

Gol 'dshleger  et  al.  {1973a)  used  a  hot  steel  ball 
with  a  diameter  of  0.4-2. 5  mm  and  smooth  surface 
finish  (<0.1  ym)  to  ignite  the  propellant.  They  noted 
that  ignition  does  not  occur  when  the  particle 
temperature  is  below  a  critical  value,  and  that 
results  are  dependent  upon  contact  between  the  hot 
particle  and  the  propellant.  They  also  proposed  a 
theoretical  model  and  obtained  numerical  solutions. 

Wise  et  al.  (1980  and  1981)  and  Law  and  Rocchio 
(1981)  conducted  similar  experiments  to  characterize 
LOVA  propellants.  They  used  hot  steel  balls  (with 
diameters  between  3.2  and  11.1  mm),  and  a  cluster  of 
propellant  slabs.  The  ignition  boundary  for 
propellants  was  determined  in  terms  of  minimum 
temperature  for  a  given  size  steel  ball,  using 
go/no-go  tests.  Their  HFCI  experiments  indicated  that 
the  binder  composition  may  be  one  of  the  most 
important  factors  in  the  determination  of  LOVA 
propellant  vulnerability. 

Kirshenbaum  et  al.  (1983)  studied  the  sensitivity 
properties  of  various  candidate  LOVA  propellants  using 
various  tests,  including  HFCI.  Their  study  showed 
that  LOVA  candidates  are  significantly  less 
susceptible  to  thermal  ignition  than  conventional 
propellants.  Caveny  et  al.  (1973)  also  investigated 
the  effect  of  additives  on  the  flammability  limits  of 
propellants,  investigated  the  effect  of  additives  on 
flammability  limits  of  propellants,  and  evaluated 
three  types  of  additives  (including  coolants,  char 
formers,  and  flame  inhibitors)  for  reducing 
flammability. 

Prior  to  the  present  study,  several  attempts  have 
been  made  to  model  the  ignition  process  encountered  in 
HFCI  experiments  (Gol'dshleger  et  al.,  1973a;  Vilyunov 
and  Kolchin,  1966;  Anderson  et  al.,  1972;  Gol'dshleger 
et  al.,  1973b;  Linan  and  Kindelan,  1981;  and  Tyler  and 
Jones,  1981).  However,  because  most  of  the  existing 
models  of  conductive  ignition  employ  numerous 
simplifying  assumptions  and  asymptotic  solution 
procedures,  their  validity  and  usage  are  limited. 
Grossman  and  Rele  (1974)  obtained  numerical  solutions 
for  the  ignition  of  cellulose  by  impingement  of  hot, 
high-velocity  spherical  particles.  They  concluded 
that  the  total  energy  of  the  metal  particle  (both 


thermal  and  kinetic)  is  important  in  determining 
whether  or  not  ignition  occurs. 

None  of  the  theoretical  models  to  date  have 
attempted  to  describe  the  foam  and/or  melt  layer  which 
has  been  experimentally  observed  on  the  propellant 
surface  in  contact  with  hot  particles.  Nor  has  the 
effect  of  binder  composition  been  studied  in  previous 
models.  A  theoretical  model  for  simulating  the  igni¬ 
tion  processes  of  solid  propellants  in  HFCI  tests  is 
needed  in  order  to  evaluate  new  propellant  formula¬ 
tions  with  respect  to  their  ability  to  resist  conduc¬ 
tive  ignition  and  thereby  reduce  system  vulnerability. 

Specific  objectives  of  this  study  are:  1)  to 
formulate  a  comprehensive  model  for  simulating  the 
ignition  of  solid  propellants  under  HFCI  situations; 

2)  to  identify  chemical  and  physical  input  parameters 
required  for  the  model;  3)  to  design  and  conduct 
simplified  melting  and  evaporation  experiments  to 
partially  validate  the  theoretical  model  under  special 
limiting  conditions;  and  (4)  to  demonstrate  the 
capability  of  the  theoretical  model  to  simulate  ice 
melting  and  evaporation  processes. 

METHOD  OF  APPROACH 


Selection  of  Particle  Geometry  and  Description  of 


In  actual  HFCI  processes,  the  hot  fragments 
(spall  particles)  generated  from  the  penetration  of 
shaped  charge  jets  of  armor  plates  are  In  different 
shapes  and  sizes.  In  order  to  understand  the  physical 
and  chemical  processes,  a  particular  particle  shape 
has  to  be  selected  in  the  formulation  of  the 
theoretical  model.  Since  no  particular  shape  can 
truly  represent  the  actual  geometry  of  spall 
fragments,  any  particle  geometry  (e.g.  spherical, 
cylindrical,  cubic,  parallel  piped)  could  be  chosen. 

The  cubic  or  parallel  piped  particles  require 
three-dimensional  and  unsteady  solution,  which  is  more 
complex  than  the  cylindrical  and  spherical  cases.  The 
spherical  particle  geometry  seems  simpler  than 
cylindrical;  however,  the  heat  transfer  process 
involved  does  not  satisfy  the  point  symmetry  condition. 
Therefore,  the  cylindrical  particle  geometry  has  been 
selected.  One  additional  advantage  of  the  cylindrical 
spall  particle  is  due  to  its  flat-bottomed  surface, 
which  can  easily  maintain  a  good  contact  with  a  flat 
propellant  surface. 

Consider  a  typical  HFCI  experiment  using  a  hot 
spall  particle  with  cylindrical  geometry.  At  time 
equals  zero,  the  spall  particle  is  placed  on  a  cold 
propellant  sample.  In  the  early  phase  of  the  process 
(Time  Period  I),  heat  is  conducted  from  the  hot 
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particle  to  the  propellant  without  any  phase  change  or 
pyrolysis  (see  Fig.  1).  As  time  progresses,  the 
temperature  of  the  propellant  increases  and  that  of 
the  particle  decreases.  Following  a  period  of  inert 
heating,  the  propellant  starts  to  decompose,  melt, 
and/or  gasify  (Time  Period  II).  Since  the  density  of 
the  metal  particle  is  much  higher  than  the  density  of 
the  decomposed  propellant,  it  displaces  the  decomposed 
propellant  and  becomes  imbedded  in  (sinks  into)  the 
propellant,  as  shown  in  Fig.  1.  The  amount  of 
imbedding  depends  on  the  temperature  and  size  of  the 
hot  particle  and  the  composition  of  the  propellant. 

The  decomposed  species  can  further  react 
exothermically  in  the  gas  and/or  condensed  (solid  or 
liquid  (foam))  phase  to  cause  ignition.  Self- 
sustained  ignition  will  occur  only  if  the  heat 
generated  by  the  exothermic  reactions  exceeds  the  heat 
losses.  The  entire  process  is  strongly  dependent  upon 
the  energy  content  of  the  hot  particle  and 
physicochemical  properties  of  the  propellant. 

Governing  Equations  and  Boundary  Conditions.  The 
mathematical  model  consists  of  governing  equations  and 
their  associated  initial  and  boundary  conditions  for 
the  hot  particle,  the  propellant,  and  the  melt  (or 
foam)  layer.  Figure  2  schematically  shows  these 
regions  and  the  coordinate  system.  The  energy 
equations  for  the  hot  spall  particle  and  the 
propellant  are  as  follows. 

Hot  spall  particle: 


k  (PsVs^  ’  r  Ir  (^s  ''  *  fe  (*s 

Propellant: 


dT. 


dT. 


aT  ^Ppr^pr^p^  ”  '^s  dz  ^^pr^pr^p^  “  r  dr  ^^pr"^  ^ 


dJ. 


*  h  <»pr  ‘  <  (2) 

During  the  HFCI  tests,  a  foam  layer  could  be 
generated,  and  the  gases  pyrolyzed  from  liquid  melt 
may  undergo  exothermic  chemical  reactions  to  generate 
heat  for  further  reactions.  In  order  to  take  the 
above  processes  into  consideration,  the  conservation 
equations  for  the  foam  region  must  be  formulated.  The 
following  assumptions  are  adopted: 

1)  The  liquefaction  process  at  the  interface  of 
the  propellant  and  liquid  melt  involves  no 
gasification; 


2)  The  gaseous  mixtures  are  produced  from  the 
base  of  the  spall  particle  or  at  the 
liquid-gas  interfaces  according  to  the 
Arrhenius  rate  law;  and 

3)  The  energy  release  process  in  the  gas  phase 
can  be  represented  by  a  single-step  forward 
reaction 

R  -  P  (3) 

where  R  represents  the  gaseous  reactant  mixtures,  and 
P  represents  the  gaseous  products. 

The  mass  continuity  equation  for  the  reactant 
species  can  be  derived  from  consideration  of  mass-flux 
balance  in  region  1  of  the  foam  layer  as 


+  m^  +  m^  -  m. 

The  mass  continuity  equation  for  the  product  species 
can  be  written  as 


K  -  ’'RX’gi  • 

-  ’'RjPg/ir*  *  ”1, 


The  source  terms  m5  and  can  be  expressed  as 

follows  according  to  the  Arrhenius  law. 


where  A.g  represents  the  Arrhenius  coefficient  of 
the  reaction  at  the  interface  of  liquid  and  gas  phases, 


The  mass  continuity  equation  for  the  liquid  phase  of 
the  foam  layer  in  region  1  can  be  given  as 


m  -  ») 

•  •  2 
9netj  38  '^P''  *  ' 

Rearranging  equations  (4)  and  (5),  one  can  obtain  the 
following  equation. 

DYo  , 

'^g.  ms  1 

By  considering  the  energy  fluxes  shown  in  Fig. 
3a,  one  can  derive  the  following  equation  describing 
the  rate  of  accumulation  of  the  gas-phase  energy. 

Bt  *  “m’I  • 

*  "g„,t 


«  o 


*  \^^P.R^b  *  ^“f.R^  ■  '^gas-liq^ij 


For  the  liquid  phase,  the  balance  of  energy  fluxes  are 
given  in  Fig.  3b,  and  the  energy  equation  can  be 


written  as 


K  [(1  -  L„(C„T;^  .  4H,;,)] 

■  -(‘"t.l  •  -  ’1 


*  ^"'^m^Ppr*'*s  ^^'Pl^melt  *  ^f,l^  ”  ^liq-so1id’***s 


■  ""gg^^P.P^b  *  ^“f.R^  *  ^gas-liq^ij 

(11) 

where 

^«f.g  -V»f.R 

(12) 

^vg  “  ^rS.R  *  "  ^R^^-v.P 

(13) 

^Pg  *  ^R^P,R  *  (1  *  ^rJ^P.P 

(14) 

both  qj^g  and  di^n-solid  positive  for  downward 
heat  fluxes.  ^ 


Initial  and  Boundary  Conditions.  In  order  to 
solve  the  temperature  distributions  inside  the  spall 
particle  and  solid  propellant,  a  set  of  Initial  and 
boundary  conditions  must  be  specified  for  equations 
(1)  and  (2).  There  are  two  sets  of  conditions:  one 
corresponds  to  Time  Period  I  (t  <  t^eit);  the  other 
corresponds  to  Time  Period  II  (t  i  tmelt)* 
initial  conditions  for  Time  Period  I  are 


T,(0)  .  T,, 

(15) 

TpO)  -  Tp, 

(16) 

The  initial  conditions  for  Time  Period  II  are  the 
temperature  profiles  solved  at  the  end  of  Time  Period 
I. 
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The  boundary  conditions  for  spall  particles  1n 
Time  Period  I  can  be  specified  readily;  some  of  the 
boundary  conditions  are  special  cases  of  those  for 
Time  Period  II.  Therefore,  the  are  omitted  here.  The 
boundary  conditions  for  the  spall  particle  of  Case  (a) 
In  Fig.  1  during  Time  Period  II  are  given  as 

“  '  ■  '•s  =  -‘s  ^  ■  "cltTs  -  T.)  *  (17) 


at  z  -  0  : 

*s  \  <Vb  *  “f.ii 

*  *8^Jas*  0  -  ’bt’il,  * 

(18) 

at  r  •  0  : 

«T  n 

TT  ■  “ 

(19) 

at  r  . 

r  r  "eff-  (t.*)CT,(t.r,.t)  -  Vt)] 

IF-  r’  ' 

• 

-  ^  (for  0  <  z  i  hj) 

(20a) 

:T7''eff2  (‘.')[T,(t.V7)  -Tgl 

®  (for  <  z  S  Lj) 

(20b) 

The  boundary  conditions  for  LOVA  propellant  In  Time 
period  II  for  Case  (a)  can  be  written  as  follows 

at  z  -  -  • 

•  ^p  ■  ^p1 

(21) 

2  •  -Lnj;  Tp  -  T^elt 

i  '"m 

(22) 

at  r  - 

Tp  -  Tp^ 

(23) 

V  ■  Vn  V  »  ») 


at  r  -  r, 
h 


tn 
<  z 


pr  dr  ^ 


radp  *  ^2(  V-  0) 

(25) 


“  ^melt 

ai. 


(if  >  0) 


at  z 


’  ’Vr  IT  "  ^3  ^"^ps  “  ■'■J  ■  ^Jad, 


r  >  r 


m 


(26) 


The  boundary  conditions  for  the  spall  particle  and 
LOVA  propellant  in  Time  Period  II  for  Cases  (b),  (c), 
and  (d)  can  be  written  in  a  similar  manner  and  are 
omitted  here  to  save  space. 


Simplifications  of  the  MathOTatic  Model. 

Although  the  above  set  of  governing  equations  can  be 
solved  with  their  boundary  and  initial  conditions, 
together  with  a  set  of  adequate  heat  transfer 
correlations,  the  calculation  of  the  above  POE  model 
is  extremely  time-consuming.  In  order  to  bypass  this 
difficulty  in  solving  many  coupled  partial 
differential  equations  (POEs)  and  intricate  boundary 
conditions,  the  theoretical  model  has  been  recast  into 
a  set  of  ordinary  differential  equations  (DDEs)  using 
the  Goodman's  (1964)  integral  method. 

In  general,  the  instantaneous  temperature 
distribution  inside  the  spall  particle  with  a 
cylindrical  shape  is  a  function  of  radial  and  axial 
positions.  During  the  interval  of  contact  with  the 
LOVA  propellant,  the  heat  fluxes  across  the  boundaries 
of  the  spall  particle  are  high  enough  to  generate 
non-negligible  temperature  gradients  inside  the 
particle.  Therefore,  lumped-parameter  analysis  cannot 
be  used  for  this  situation.  To  obtain  approximate 
solution  of  the  HFCI  model  based  upon  POE  formulation, 
the  temperature  profile  is  assumed  to  obey  the 
following  polynomial  form  (see  Fig.  4). 


T-(t,r,z) 

Tu(t) 


[1  *  Cj(t)r^  *  C2(t)r^]  • 


[1  *  CjdXz  -  zj  *  C^(t)(2  - 


*  C5(t)(2  -  ♦  Cj(t)(z  -  z„)^] 


(27) 


where  Ty{t)  represents  either  the  temperature  of  the 
uniform-temperature  zone  in  the  spall  or  the  maximum 
instantaneous  temperature  when  the  uni  form- temperature 
zone  shrinks  to  zerP'.  The  height  z^,  represents  either 
that  of  the  midpoint  of  the  uniform-temperature  zone 
or  that  of  the  maximum  temperature  point.  can  be 
calculated  from  the  following  equations 


Ir.  ^  ,  .1  for  finite  uniform-, 

?^®thg  ^^s  ”  ®thj^^  temperature  zone 


after  uniform-temperature 
zone  disappears 


(28)b 


Ci{t)  through  C5(t)  are  time-dependent  coefficients  to 
be  determined  by  various  boundary  and  smooth 
conditions  discussed  later. 

Following  Goodman's  (1964)  integral  method,  by 
integrating  equation  (1)  twice  with  respect  to  z  and  r 
and  applying  the  polynomial-  form  of  equation  (27),  the 
POE  for  the  spall  particle  is  reduced  to 


where  the  parameters  P(t)  and  S(t)  are  defined  as 


P(t)  ■  [ej(t)  ♦  e^(t)  ♦  e^(t)  ♦  9„(t)]  (30) 

S(t)  =  i  ^  [r/ -  (r,  -  *  I  Cj- 

*  *  2(2Ci  ‘  ■ 


r  cc-s  -  *  r  [f'-s  - 


•  r  cc-s  -  -  VI  *  r  [c-s  -  -  V^l 


In  equation  (30),  e,.  Or.  and  8,.  represent  the 
thermal  energy  content  in  the  bottom,  top,  and  lateral 
thermal  penetrated  zone  of  the  spall  particle;  their 
mathematical  definitions  are  as  follows: 


(32) 


es(‘.-r  I 

-^0 


S  ''th^(t) 


C^p^T^2irrdrdz 


Cjp5T^2irrdrd2  (33) 


e,  (t) 


f  *  '•"T  f  =* 


s 

^-«th  (‘)  <”) 


The  thermal  energy  storage  in  the  uniform  temperature 
zone,  6u(t),  can  be  expressed  as 


•-s-^th  % 

I  s  tny  |.  s 


The  above  four  0s  are  expressed  In  terms  of  thermal 
wave  penetration  depths  and  constants  Cx  through  Cg. 
These  six  constants  are  solved  from  a  set  of  coupled 
linear  algebraic  equations  obtained  from  the  boundary 
and  smooth  conditions.  Their  explicit  forms  are  given 
by  Kuo  et  a1.  (1985). 

To  evaluate  the  thicknesses  of  the  thermal  waves 
near  the  top  and  bottom  of  the  spall  particle,  the 
following  equations,  deduced  from  the  expression 
obtained  by  Goodman  (1964),  can  be  used. 


(t) 


/55  (p 


eff 


1C 


’mr. 


•tc 


-  Tg)dt} 


1/2 


(36)a 


6th.(t)  -  /So  [- 


r 


O^pall  » 


^spall 


(t)dt] 


1/2 


(36)b 


where  the  subscript  TC  stands  for  top  center  location 
of  the  spall  particle.  The  thickness  of  the  thermal 
wave  penetration  depth  in  the  lateral  direction  can 
also  be  given  in  a  similar  form. 
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Using  a  third-degree  polynomial  in  z*  to  repre¬ 
sent  the  temperature  profile  in  the  solid  propellant 
and  integrate  equation  (2)  with  respect  to  z*.  the 
time  variation  of  the  thermal  wave  penetration  depth 
(the  thickness  of  the  heated  zone)»  6pr>.  can  be 
determined  from 


u> 

-f - 6  V  •  V'‘> 

The  initial  condition  for  6pr  is 

fip^(O)  -  0  (38) 

After  6pr(t)  is  solved,  qp^(t)  can  be  calculated 
from  ^ 

<<*>  ■  '3*) 

Heat  Transfer  Correlations  Employed  in  the 
Theoretical  ModeK  To  complete  the  model  formulation, 
the  following  heat -transfer  correlations  are  presented 
as  a  part  of  physical  input  to  the  model. 

a.  Heat  Transfer  at  the  Base  of  Spall  Particle. 
The  rates  of  heat  transfer  from  the  spall  particle 
base  to  th^  liquid  and  gas  bubbles  in  the  foam  layer 
(q^^q  and  qg^j)  can  be  evaluated  from 

’ilq  ■ 


’Jas  ■  "base  l^B  '  («) 

When  the  average  temperature  at  the  base  of  the 
spall  particle  (Tsb)  is  higher  than  the  saturation 
temperature  of  the  liquid  melt,  the  boiling  heat 
transfer  correlation  should  be  used.  According  to  the 
empirical  correlations  of  Jakob  and  Hawkins  (1957)  and 
Holman  (1977),  hpase  related  to  the  temperature_ 
difference  between  150  and  saturation  temperature  Tsat 
by  the  following  equation. 


(42) 


,  ‘  >5770. 

^'base  - 

,5-56  (Ts8  -  T„t)  ‘5”"  ^  2-365  «  10 

This  equation  was  developed  from  boiling  heat-transfer 
data  of  water  at  one  atmospheric  pressure  on  a 
horizontal  plate.  Strictly  speaking,  this  correlation 
is  not  completely  suitable  for  the  downward  heating 
from  the  spall  particle  base  to  the  foam  layer,  but  in 
the  absence  of  better  correlations,  it  is  temporarily 
adopted  and  used  with  caution.  According  to 
Cheung  (1985),  the  boiling  heat  transfer  for  the 
downward-heating  case  may  follow  the  same 
temperature-difference  dependence,  but  the  constants 
could  be  approximately  two  to  three  times  higher  than 
the  upward  heating.  In  view  of  the  confined  space 
available  between  the  base  of  the  spall  particle  and 
the  solid  propellant,  the  gas  bubble  may  be  trapped 
momentarily  in  this  foam  layer  region.  If  this 
happens,  the  material  in  the  foam  layer  could 
accumulate  thermal  energy  and  reduce  the  rate  of 
boiling  heat  transfer.  Because  the  downward-heating 
effect  may  compensate  for  the  confined-space  effect, 
equation  (42)  is  used  at  present  for  the  simulation  of 
melting  and  evaporation  of  ice  to  be  discussed  in  a 
later  section. 


After  the  base  temperature  of  the  spall  particle 
has  dropped  below  the  boiling  temperature,  gas  bubbles 
could  still  be  generated  from  the  spall  surface  (due 


to  transient  heating  of  the  foam  layer)  until  dspan 
decreases  to  the  level  of  q,JJ^n  .  q,l^„  can  be  ^ 
calculated  from  the  following  equation  based  upon  the 


studies  of  Chang  (1959),  Zuber  (1959),  and  Lienhard 


and  Dhir  (1980). 


0-177  Pgh,g 


1/4 


(43) 


where  g  is  the  gravitational  acceleration. 

The  average  velocity  of  bubbles  leaving  the  heated 
surface  area  is 


,  1/4 

Vjj  -  0.59  [gOj  (pj  -  Pg)/Pt  ]  (44) 

which  is  based  upon  the  relationship  reported  by 
Rohsenow  and  Choi  (1969).  The  average  bubble  diameter 


at  detachment  from  the  heating  surface,  based  upon  the 
empirical  formula  of  Fritz  (1935),  is 


Dp  - 


5157^ 


(45) 


where  B  is  approximately  40®  and  Cj  -  0.0148  for  water 
in  the  ice  melting  and  evaporation  simulation. 

The  void  fraction  at  the  base  of  the  spall 
particle,  ifB,  can  be  evaluated  from  the  following 
approximation. 


(46) 


where  a  is  a  stretching  constant  in  the  order  of  one, 
specific  heat-flux  ratio  at  which  versus 
has  a  point  of  inflection.  This  equation 
ed  based  upon  the  fact  that  it 
satisfies  the  limiting  case  of  film  boiling,  which 
requires  that  'Vq  *  1  when  qjpgil  »  It  also 

When  Tsb  is  lower  than  the  saturation  tempera¬ 
ture  of  the  liquid  melt,  the  heat-transfer  coefficient 
between  the  fluid  in  the  foam  layer  and  the  bottom 
surface  of  the  spall  particle  can  be  evaluated  from 
the  following  Kercher's  empirical  correlation  taken 
from  the  General  Electric  Heat  Transfer  Data  Book 
(1979). 

r*  A  1  n  u  n  I 

^mixture  ,  m^ 

W - (%) 

where  A^iixture  calculated  from 


and  b  i§  a 

^spall'^^tnin 
IS  formulat 


^mixture  "  '^'s^vapor  *  "  *^8^  Miq 


(48) 


(j»l  and  ^2  equation  (47)  are  correction  factors 
which  depend  upon  the  ratio  of  bubble  separation 
distance  (Xp),  bubble  diameter,  and  the  Reynolds 
number  (■PgVbDb/i^g)  •  above  equation  was  initially 
developed  for  gases  to  flow  through  an  array  of 
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circular  holes  of  a  porous  plate,  located  at  the  top 
of  a  horizontal  gap.  Numerous  jets  impinge  on  a 
bottom  horizontal  plate  located  a  small  distance  from 
the  top  porous  plate.  The  bubbles  generated  at  the 
base  of  the  spall  can  be  regarded  as  gases  injected 
downward  through  a  horizontal  porous  plate.  Even 
though  the  locations  of  bubbles  vary  from  time  to 
time,  as  long  as  the  bubble  sizes  and  separation 
distances  are  properly  calculated,  the  flow  agitation 
and  heat  transfer  processes  are  quite  similar  to  those 
of  a  porous  plate.  The  average  separation  distance 
between  bubbles  is  calculated  from  the  bubble  size  and 
the  void  fraction  in  the  foam  layer  by  the  following 
equation. 


X 


n 


(49) 


The  heat-transfer  coefficient  given  in  equation  (47) 
is  used  to  determine  the  energy  fluxes,  qg^j  and  q'-j-jq. 
until  Tsb  decreases  to  T^elt*  ^ 


b.  Heat  Transfer  at  the  Interface  Between  Liquid 
Melt  and  3o1id  Material.  The  heat-transfer  process  in 
the  foam  layer  at  the  interface  with  solid  material  is 
governed  by  the  condensation  of  bubbles  at  the 
interface.  Levenspiel's  correlation  of  steam-bubble 
condensation  (Soo,  1967)  has  been  adopted  to  calculate 
the  instantaneous  heat-transfer  coefficient  and  the 
heat  flux  qiia-soiid*  transfer  coefficient 

can  be  written  as 


h^(t)  -  3.54  (50) 

To  adequately  calculate  the  values  of  hc(t),  the 
effective  bubble  diameter  is  used  in  the  above 
equation  to  replace  the  instantaneous  diameter  of  a 
single  bubble.  During  the  condensation  process,  the 
rate  of  decrease  of  bubble  diameter  can  be  described 
by  the  following  equation  (Soo,  1967). 

d  tn  D.(t) 

- 5j5_.7.08  [yt.O)  -Tg^(t)]  (51) 


Both  the  propellant  or  ice  surface  temperature  and  the 
bulk  temperature  of  the  gas  bubble  are  functions  of 
time.  However,  experimental  evidence  indicates  that 
the  bubble  lifetime  is  much  shorter  than  the 
foam-layer  action  time;  hence,  Tg  can  be 
treated  as  a  constant  in  estimating  single  bubble 
lifetime.  Setting  the  propellant  surface  temperature 
equal  to  Tmaif,  equation  (51)  can  be  integrated  with 
respect  to  time  to  give 

Djj(t)  -  Djj(O)  exp[-7.08  (T^^-  t]  (52) 

Using  the  same  functional  form,  the  effective  bubble 
diameter  can  be  expressed  as 

0.  (t)  -  D.  exp[-kt]  (53) 

°eff  ° 

where  O5  (the  average  diameter  at  detachment)  is 
calculated  from  equation  (45)  and  k~l  represents  the 
characteristic  decay  time.  The  physical  meaning  of 
the  effective  bubble  diameter  is  associated  with  the 
averaged  void  fraction  of  the  surface  area  due  to 
bubble  impingement  on  the  propellant  surface.  As  the 
spall  particle  sinks  into  the  ice  block,  the  particle 
is  cooled  by  transferring  the  heat  to  its  adjacent 
material.  The  boiling  process  becomes  less  intense 
and  generates  fewer  gas  bubbles  to  impinge  on  the  ice 
block.  Therefore,  the  effective  diameter  decreases 
with  respect  to  time.  The  k  value  used  in  equation 
(53)  is  in  the  order  of  0.2  sec'^.  To  account  for  the 
above  effect  caused  by  the  reduction  of  bubble  number 
density  in  terms  of  effective  bubble  diameter,  the 
decay  time  (k"l)  used  in  the  computation  is 
considerably  longer  than  the  bubble  lifetime. 

c.  Heat  Transfer  at  the  Top  and  Lateral  Surfaces 
of  the  Spall  Particle,  the  top  surface  of  a  spall 
particle  can  be  considered  as  a  horizontal  plate  where 
the  heat-transfer  mechanism  is  due  mainly  to  natural 
convection  (McAdams,  1954;  and  Goldstein  et  al.,  1973). 
The  Nusselt  number  is  expressed  as  a  function  of  the 
Rayleigh  number  in  the  commonly  used  McAdams'  (1954) 
correlation. 

The  lateral  surface  of  the  spall  particle  in 
contact  with  the  foam  layer  experiences  natural 
convection  or  conduction  since  the  velocity  of  the 
melt  layer  is  quite  small.  On  the  other  hand,  the 
lateral  surface  of  the  spall  particle  which  is  exposed 
to  the  air  experiences  forced  convection  due  to  the 
flow  of  pyrolyzing  gaseous  species  above  the  melt 


layer.  In  order  to  consider  these  two  different 
mechanisms,  it  is  necessary  to  determine  the  different 
heat-transfer  coefficients  caused  by  forced  and  free 
convection. 

For  laminar  flows,  the  conventional  correlation 
for  parallel  external  flow,  along  the  axis  of  a 
cylinder,  is  used  to  relate  the  local  Nusselt  number 
to  the  local  Reynolds  and  Prandtl  numbers  (Holman, 
1972).  For  turbulent  flows,  Schlichting's  (1968) 
local  friction  coefficient  is  used  to  relate  Reynolds 
numbers  with  the  local  Nusselt  number. 

Application  of  the  HFCI  Model  for  Simulating  the 
Melting  and  Evaporation  of  Ice  upon  Direct  contact 
with  a  Hot  Elemenf 

In  order  to  verify  the  theoretical  model 
developed  for  predicting  HFCI  phenomena,  a  set  of 
input  parameters  and  functions  is  required.  At 
present,  although  some  measurements  are  being 
conducted  by  Miller  (1986),  many  critical  parameters 
for  LOVA  propellants  are  not  fully  characterized.  In 
view  of  this  situation,  experiments  must  be  conducted 
in  a  simpler  system  to  simulate  HFCI  phenomena. 

One  obvious  way  to  simulate  the  conductive 
heating,  melting,  and  evaporating  processes  of  HFCI 
tests  is  to  heat  an  ice  block  with  a  hot  metal  rod. 

In  such  an  experiment,  uncertainties  in  chemical 
kinetics  in  the  melt  layer  are  avoided,  and  data  from 
experiments  can  be  used  to  validate  the  theoretical 
model  under  the  limiting  condition  without  involving 
chemical  reactions.  The  reason  ice  was  selected  for 
melting  and  evaporation  simulation  is  due  to  the 
well-known  boiling  heat  transfer  characteristics  of 
water,  as  well  as  its  thermal  properties.  The  ice 
cube  is  also  readily  available  for  experiments. 

The  regression  velocity  of  the  ice  surface,  v^, 
is  computed  from  the  heat-flux  balance  at  the 
interface  between  the  foam  layer  and  ice  by  the 
following  equation. 


Major  and  minor  unknowns,  together  with  the 
numbers  of  equations  for  their  determination,  are 
listed  in  Tables  1  and  2. 
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Table  1.  Major 

Unknowns 

and  Equations 

Region 

Unknown 

Equation 

Spall 

TS 

(27) 

Particle 

P 

(29) 

Tu 

(30). (31-35) 

Foam 

Yr 

-0 

Layer 

(10) 

(11) 

Vm 

(54) 

4' 

(4) 

Ice 

•Spr 

(37) 

Table  2.  Minor  Unknowns  and  Equations 


Region 

Unknown 

Equation 

Unknown 

Equation 

00 

(32) 

0L 

(34) 

Spall 

Particle 

0T 

(33) 

©u 

(35) 

<5thT 

(36)a 

«thB 

(36)b 

^thL 

(36)a  4  (36)b 

2m 

(28)a  or  (28)b 

AHf.g 

(12) 

Db 

(45) 

Foam 

Layer 

Cvg 

(13) 

Cpg 

(14) 

^’base 

(42)  or  (47) 

qiiq 

(40) 

he 

(50) 

^5as 

(41) 

'J'B 

(46) 

Ice 

Qpr 

(39) 

TEST  SETUP  AND  RESULTS  FOR  ICE  MELTING  AND  EVAPORATION 
EXPERIMENTS 

An  ice-melting  and  evaporating  experiment  has 
been  devised  and  conducted.  To  be  compatible  with  the 
theoretical  model,  a  steel  cylindrical  rod  with  a 
diameter  of  1.27  cm  and  a  length  of  2.54  cm  is 
connected  at  one  end  to  a  10.16  cm-long  guiding  steel 
rod  with  a  diameter  of  0.3175  cm.  Four  Copper- 
Constantan  thermocouples  with  bead  sizes  of  0.025  cm 
and  wire  diameters  of  0.0075  cm  were  soldered  to 
different  locations  of  the  heating  element  (shown  in 


Fig.  5).  In  these  tests,  the  heating  element  was 
heated  inside  an  oven  with  preselected  temperature. 

The  heating  element  was  placed  on  top  of  an  ice  cube, 
and  the  temperature-time  variations  at  different 
locations  were  recorded.  A  typical  set  of  recorded 
temperature-time  traces  is  shown  in  Fig.  5. 

In  addition  to  temperature-time  measurements,  a 
video  movie  camera  was  used  to  record  the  ice  melting 
and  evaporating  phenomena.  The  instantaneous  position 
of  the  heating  element  was  recorded  on  the  video 
camera.  The  displacement  and  velocity  of  the  heating 
element  were  determined  by  reading  the  instantaneous 
top  surface  location  of  the  heating  element  against  a 
stationary  scale  and  the  digital  time  counter  on  the 
video  screen.  Time  variations  of  the  measured  velo¬ 
city  of  the  heating  element  are  compared  with  theoret¬ 
ical  results  in  a  later  section. 

In  order  to  observe  the  flow  pattern  of  the  gas 
and  steam  surrounding  the  heating  element,  a  Schlieren 
optical  system  {Settles,  1970;  and  Kuo  et  al.,  1985) 
was  used.  A  typical  color  Schlieren  photograph 
obtained  from  an  ice-melting  and  evaporating 
experiment  is  shown  in  Fig.  6.  Several  interesting 
observations  are  made  from  the  video  and  temperature 
records  of  these  experiments: 

1)  Although  all  four  temperature-time  traces  in  Fig. 

5  showed  similar  variations,  there  are  some 
differences  among  these  traces.  The  temperatures 
on  the  main  body  of  the  heating  element  decay  at 
a  slower  rate  than  those  on  the  guiding  rod. 

This  is  due  to  the  fact  that  the  surface  area  per 
unit  volume  is  much  larger  on  the  guiding  rod 
than  that  of  the  main  body. 

2)  Based  upon  the  slope  changes  of  the  temperature 
decay,  the  overall  process  can  be  divided  into 
the  three  time  zones  shown  in  Fig.  5.  In  the 
initial  time  zone,  the  temperature  of  the  heating 
element  decays  slowly  with  time  since  the  element 
was  taken  from  the  oven  and  cooled  by  natural 
convection  in  air  and  radiation  to  the 
surrounding.  A  small  amount  of  energy  is  lost 
from  the  guiding  rod  to  a  room-temperature  clamp 
by  conduction.  In  the  second  time  zone,  a 
significant  amount  of  steam  was  generated  at  the 
contact  zone  of  the  heating  element  and  the  ice. 
The  flow  around  the  cylindrical  rod  was  highly 
turbulent.  Associated  with  steam  generation,  the 
temperatures  at  all  thermocouple  locations 
decrease  drastically  in  1.2  seconds,  as  shown  in 
Fig.  5.  In  the  last  time  zone,  steam  generation 
is  gradually  replaced  by  ice  melting  and  slow 
sinking  of  the  heating  element  into  the  ice  block. 
Eventually,  the  heating  element  stopped  at  a 
terminal  position  inside  the  ice  block. 
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3)  The  downward  velocity  of  the  heating  element  was 
highest  upon  initial  contact  with  the  ice  block, 
and  then  decayed  monotonical ly  to  zero  velocity. 

4)  From  the  color  Schlieren  photographs,  one  can  see 
clearly  that  the  steam-jet  velocity  is  mostly  in 
the  radial  direction  at  the  initial  time  interval. 
The  steam  jet  then  changed  from  radial  outward 
direction  to  vertical  upward  direction.  Associ¬ 
ated  with  steam-jet  action,  the  flow  field  sur¬ 
rounding  the  cylinder  is  highly  turbulent. 
Following  the  decay  of  steam  generation,  the  flow 
became  laminar  at  the  later  stage. 

COMPARISON  OF  NUMERICAL  RESULTS  WITH  TEST  DATA  OF  ICE 
MELTING  AND  EVAPORATION 

To  partially  verify  the  HFCI  model,  numerical 
solutions  were  compared  with  experimental  data 
obtained  from  an  ice-melting  and  evaporating 
experiment.  Input  data  for  this  numerical  simulation 
is  given  in  Table  3.  Figure  7  presents  a  comparison 
of  the  calculated  sinking  velocity  of  the  spall 
particle  with  the  measured  data.  It  is  quite  obvious 
that  agreement  between  calculated  results  and  measured 
data  is  reasonably  close.  The  trend  of  sharp  decline 
followed  by  a  much  slower  rate  of  sinking  is  exhibited 
in  both  theoretical  and  experimental  results.  Figure 
8  shows  the  comparison  of  calculated  spall  particle 
trajectory  with  experimental  data.  The  agreement  is 
not  extremely  close;  however,  the  calculated  trend  of 
traveling  distance  variation  with  time  and  the 
magnitude  of  hp  are  not  far  from  the  measured  data. 


Table  3. 

Input  Data  Used  in  the 

Ice-Melting  and 

Evaporation  Simulation 

Region 

Parameter 

Value  1 

Parameter 

Value 

Spall 

^s 

0.635x10-2 

6.39x10-1 

Particle 

Ls 

2.54x10-2 

^-s 

4.34x102 

Ps 

7.83x10^ 

Tsi 

4.90x102 

Pt 

9.87x102 

4.52x10-2 

Foam 

4.18x103 

A»f.l 

-1.59x107 

Layer 

^pg 

1.87x103 

^«f.g 

-1.34x107 

Cvg 

1.41x103 

Ice 

Sr 

2.04  X  103 

^pr 

2.21 

Ppr 

9.13x102 

-1.55x107 

From  the  comparison  of  the  calculated  and  measured 
temperature-time  traces  at  the  thermocouple  TC2 
location  shown  in  Fig.  9,  it  is  evident  that  the  two 
traces  are  very  close.  In  the  same  figure,  calculated 
traces  for  the  average  temperature  at  the  base  of  the 
spall  particle  '(T33]  and  the  maximum  spall  particle 
temperature  (Ty)  are  plotted.  It  is  interesting  to 
note  that  the  temperature  difference  between  Ty  and 
Tsb  is  much  larger  than  that  between  Ty  and  TC2.  This 
implies  that  a  significant  temperature  gradient  exists 
near  the  base  of  the  spall  particle  because  of  the 
high  heat  flux  leaving  the  spall  particle  base  region. 


The  temperature  gradient  in  the  radial  direction  in 
the  middle  of  the  cylindrical  spall  particle  is  nearly 
negligible  due  to  relatively  low  radial  heat  fluxes 
and  high  thermal  conductivity  of  the  material. 

Figure  10  shows  the  calculated  radial  temperature 
profiles  at  the  top  surface,  middle  section,  and 
bottom  surface  of  the  spall  particle  for  various  times. 
Most  of  the  profiles  are  quite  flat,  but  there  is  some 
gradient  in  the  radial  direction,  especially  during 
the  early  part  of  the  event.  Calculated  temperature 
profiles  in  the  ice  are  shown  in  Fig.  11.  It  is  clear 
that  the  thermal  wave  penetration  depth  increases 
monotonically  with  respect  to  time,  while  the  surface 
temperature  is  maintained  at  a  fix^  va]ue. 

The  calculated  time  variations  of  and 

are  shown  in  Fig.  12.  It  is  useful  to  point  out  that 
^soall  significantly  higher  than  qi^  for  a  long 
penoa  of  time.  This  indicates  that  a  substantial 
amount  of  energy  is  carried  away  by  evaporation  of 
liquid  melt, in  the  foam  layer.  The  sharp  declines  of 
^spall  *^pr  initial  time  period  dictate  the 

rapid  decrease  of  sinking  velocity  with  respect  to 
time. 

SUMMARY  AND  CONCLUSIONS 

(1)  A  comprehensive  theoretical  model  for  simulating 
hot  fragment  conductive  ignition  processes  of 
LOVA  propellants  has  been  formulated.  This  model 
simulates:  heat-transfer  processes  between  the 
spall  particle  and  propellant;  formation  of  the 
liquid-melt  or  foam  layer;  propellant  pyrolysis 
processes;  displacement  of  the  hot  spall 
particle;  and  propellant  ignition  or  quenching  of 
the  spall  particle. 

(2)  The  partial  differential  equations  governing  the 
transient  heat-transfer  processes  in  the  spall 
particle  and  LOVA  propellant  have  been 
successfully  converted  to  a  set  of  ordinary 
differential  equations.  This  leads  to  a  major 
reduction  in  computational  time. 

(3)  The  HFCI  model  has  been  applied  to  the  ice 
melting  and  evaporation  processes  under  the 
direct  contact  with  a  hot  element.  Experimental 
tests  were  conducted  with  temperature 
measurements  and  flow  visualization  using  a  color 
Schlieren  system. 


(4)  Good  agreement  between  calculated  results  and 
test  data  of  ice-melting  and  evaporation 
experiments  shows  that  the  HFCI  model  is  capable 
of  predicting  the  major  heat-  and  mass-transfer 
characteristics  in  the  melting  and  evaporation  of 
a  solid  with  known  thermal  properties  and  boiling 
characteristics  of  its  liquid  melt. 

(5)  During  the  initial  contact  between  the  spall 
particle  and  ice  block,  a  strong  steam  jet  in  the 
radial  direction  is  generated.  As  time 
progresses,  the  steam  jet  changes  from  the  radial 
to  the  vertical  direction  as  its  strength  decays. 
A  substantial  amount  of  energy  is  carried  away  by 
the  steam  jet  generated  from  evaporation  of 
liquid  melt  in  the  foam  layer. 

(6)  The  input  parameters  and  types  of  correlations 
required  for  studying  the  HFCI  processes  of  a 
given  solid  propellant  can  be  identified  from  the 
theoretical  model  presented  in  this  paper.  The 
full  validation  of  the  HFCI  model  cannot  be 
performed  until  a  set  of  input  parameters  and 
correlations  become  available. 
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VARIOUS  CASES  CONSIDERED 

TIME  PERIOD  I  (BEFORE  THE  FORMATION  OF  A  MELT  LAYER) 

HOT  METAL  CYLINDER 

solid'  PROPELLANT 

TIME  PERIOD  n  (AFTER  THE  FORMATION  OF^A  MELT  LAYER) 


CASE  (b)  :  h,>  hp 


CASE(c) :  h^>hp>Ls 


CASE  (d)  :  hp>h^>L5 


Fig.  1  Physical  Event  of  Hot  Fragment  Conductive 
Ignition  (HFCI)  Processes  Represented  by 
Various  Cases  in  Two  Different  Time  Periods 


I 


SPALL  PARTICLE 


I 


a 


0 


Oir7r7rT>'iyr^'L!!y  1"  ^ 


SOLID  •  •• 
PROPELLANT- 


^gas— liq 


-.-•N 


A.- 


•L 


m 


■  |^^f,g't*PpgTgi^.  pQ 


Fig.  3a  Energy  Fluxes  Associated  with  the  Gas  Phase 
In  Region  1  (The  Thickness  of  Foam  Layer  is 
Highly  Exaggerated) 
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Fig.  3b  Energy  Fluxes  Associated  with  the  Liquid 

Phase  in  Region  1  (Thickness  of  Foam  Layer  is 
Highly  Exaggerated) 
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Fig.  7  Comparison  of  Calculated  Sinking  Velocity  of 
the  Spall  Particle  with  the  Measured  Data 
Obtained  from  Ice-Melting  and  Evaporating 
Experiment 
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Comparison  of  the  Measured  Temperature-Time 
Trace  of  TC2  with  Calculated  Temperature- 
Time  Traces  at  the  Spall-Particle  Base,  TC2 
and  Maximum  Temperature  Locations  of  the  Spall 
Particle 
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Fig.  12  Calculated  Time  Variation  of  Heat  Fluxes 
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Modeling  of  Hot  Fragment  Conductive  Ignition  Processes 

of  LOVA  Propellants 
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ABSTRACT 

A  comprehensive  theoretical  model  and  an  efficient  numerical  program  have  been  developed  to 
simulate  the  hot  fragment  conductive  ignition  (HFCI)  processes  for  characterizing  the  degree  of 
vulnerability  of  various  gun  propellants.  In  the  formulation  of  the  theoretical  model,  three  sets 
of  governing  equations  and  their  boundary  conditions  were  derived  for  different  regions  consisting 
of  the  hot  spall  particle,  the  LOVA  propellant,  and  the  foam  layer  produced  by  liquefaction, 
pyrolysis,  and  decomposition  of  the  propellant.  In  terms  of  chemical  kinetic  scheme,  a  two-step 
reaction  model  la  proposed,  based  upon  DSC  experiments.  The  model  contains  an  endothermic 
decomposition  process  followed  by  an  exothermic  pyrolysis  of  the  propellant.  Calculated  results 
compared  well  with  experimental  data  in  temperature-time  traces,  trajectory  of  spall  particles 
sinking  into  the  LOVA  propellant,  and  teralnal  position  of  the  particle  for  quenched  cases. 
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Universal  gas  constant,  831  4. ii  J/kg-mole-K 
radius  of  spall  particle,  n 
boiling  temperature  of  liquid  melt,  K 
gas  temperature,  K 

bulk  temperature  of  melt  or  foam  layer,  K 
melting  temperature  of  LOVA  propellant,  K 
temperature  of  LOVA  propellant,  K 

average  surface  temperature  of  LOVA  propellant  exposed  to  air,  K 
temperature  of  spall  particle,  K 

radial  velocity  of  melt  or  foam  layer  In  Region  1 ,  m/s 

sinking  velocity  of  spall  particle,  m/s 

axial  distance  above  base  of  spall  particle,  m 

axial  distance  below  instantaneous  surface  of  propellant  under  foam  layer,  m 


Greek  Symbols 

a  thermal  dlffuslvlty,  m^/a 

6pr  thermal  penetration  depth  from  propellant  surface,  m 

X  thermal  conductivity,  W/m-K 

p  density,  kg/m3 

i|i  average  porosity  (void  fraction)  of  foam  layer  In  Region  1 

u  dynamic  viscosity,  kg/m*a 


8 

i 

P 

malt 
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• 
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Subscripts 

gas-phase 

initial  condition 

liquid-phase 

melt  liquid  melt 

propellant 

room  condition 

Region  1  of  foam  layer 


Superscript 
■  T  -  Tref 


INTRODUCTION 

In  the  selection  of  low  vulnerability  amminition  (LOVA)  propellants,  go/no-go  ignition  tests 
with  hot  metallic  elements  are  often  conducted  to  simulate  hot  fragment  conductive  Ignition  (HFCI) 
processes.  These  hot  fragments  can  oe  considered  as  those  generated  from  the  penetration  of  armor 
plates  by  shaped  charge  Jets  or  kinetic  energy  penetrators.  One  way  to  negate  the  threat  of 
propellant  ignition  by  hot  spalls  is  to  use  a  propellant  which  la  resistant  to  conductive  Ignition. 

A  number  of  Investigations  on  this  subject  have  been  conducted,  both  theoretically  and 
experimentally,  in  recent  years  [1-12],  Since  a  review  of  pertinent  work  is  given  in  the  paper  by 
Kuo  et  al.  [12],  detailed  reviews  are  omitted  here.  Only  major  findings  of  HFCI  processes  are 
listed  below. 

1.  Binders  of  nitramine  composite  propellants  were  found  to  have  a  strong  effect  on  conductive 
ignition.  Some  binders  can  act  as  fire  retardant  coolants. 

2.  The  relative  susceptibility  of  LOVA  propellants  to  ignition  by  spall  particles  can  be  determined 
from  ignition  map  based  upon  the  plot  of  the  Initial  temperature  versus  weight  of  the  spall 
fragment. 

2.  The  contact  resistant  between  hot  particle  and  propellant  slgnif Icantlyaffects  the  delay  time 
for  onset  of  runaway  ignition. 
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4.  Binders  which  endothermically  decompose  under  acid  catalysis  are  more  desirable  for  uOVA 
propellant  binder  ingredients. 

3.  The  Ignition  temperatures  of  LOVA  propellants  are  similar  to  their  nltramlne  fillers. 

6.  For  certain  LOVA  propellants,  liquefaction  and  bubble  formation  were  observed  during  conductive 
heating.  This  suggests  that  the  thermal  Insulation  properties  of  binder  products  may  be 
Important  In  the  determination  of  Ignition  sensitivities  of  LOVA  propellants. 

7.  Based  upon  surface  thermocouple  measurements  during  HFCI  tests,  the  Initial  condensed  phase 
reactions  are  endothermic.  The  endothermlcity  of  these  reactions  showed  good  correlation  with 
propellant  sensitivity  for  thermal  Ignition. 

8.  The  minimum  spall  particle  temperature  (for  a  fixed  spall  particle  weight)  depends  strongly  on 
binder  composition. 

9.  Chemical  reactions  leading  to  Ignition  of  certain  LOVA  propellants  during  HFCI  experiments  can 
be  idealized  as  a  two-step  sequence  of  global  reactions  in  which  an  endothermic  reaction  Is 
followed  by  an  exothermic  reaction. 

Several  attempts  were  made  to  model  the  hot  fragment  conductive  Ignition  processes  [9-11]; 
however,  these  models  have  serious  limitations  In  simulating  actual  HFCI  processes.  The  model 
proposed  In  Ref.  9  Is  strictly  one-dimensional,  and  allows  no  space  for  the  products  In  the  gas  or 
liquid  phases  to  exit  at  the  Interface  between  propellant  samples  and  hot  Inert  fragments.  The 
model  proposed  In  Ref.  10  requires  no  depletion  of  propellant  when  the  hot  fragment  penetrates  the 
combustible  material.  Also,  no  surface  reaction  Is  allowed  in  the  model.  An  Interesting 
theoretical  analysis  of  the  Ignition  of  reactive  solids  by  direct  contact  with  a  hot  inert  body  was 
performed  by  Llnan  and  Kindelan  [11].  Their  analysis  considered  both  cylindrical  and  spherical 
geometries.  Asymptotic  solutions  for  large  activation  energies  were  obtained.  The  solution 
Indicated  the  existence  of  an  Ignition  boundary  beyond  which  runaway  Ignition  occurred.  Although 
Che  theoretical  results  Indicated  the  correct  trend  In  terms  of  Inert  particle  initial  temperature 
and  particle  size,  their  solutions  were  not  compared  with  any  experimental  data.  Furthermore,  the 
assumption  of  high  activation  energies  are  Invalid  for  various  LOVA  propellants.  There  was  no 
provision  for  the  development  of  a  foam  layer  at  the  Interface  region  between  hot  fragment  and 
propellant. 

A  comprehensive  theoretical  model  was  recently  developed  by  Kuo  et  al.  [12]  to  describe  the 
heat  transfer  processes  In  hot  fragment  and  solid  propellant,  as  well  as  the  development  of  a  foam 
and/or  melt  layer  which  has  been  experimentally  observed  at  the  interface  zone  between  solid 
propellant  surface  and  hot  spall  particle.  In  this  model,  the  liquefaction,  surface  reaction,  and 
gas-phase  reactions  In  the  foam  layer  were  considered.  The  effect  of  binder  composition  can  also  be 
studied  by  using  measured  kinetic  data  from  separate  OSC  experiments  for  LOVA  propellants  with 
various  binders.  Thus,  the  ability  of  resisting  conductive  ignition  of  new  propellant  formulations 
can  be  evaluated  for  achieving  reduced  system  vulnerability. 

Although  the  model  presented  In  Ref.  12  Includes  both  reacting  and  non-reacting  cases,  only  the 
non-reacting  part  of  the  model  was  verified  by  Ice-meltlng  and  evaporation  experiments.  The  main 
objective  of  this  paper  la  to  verify  the  chemically  reacting  part  of  the  model.  The  latest  reaction 
mechanism  and  data  obtained  by  Killer  et  al.  [8]  have  been  Incorporated  Into  the  model.  Numerical 
simulation  of  HFCI  processes  of  certain  LOVA  propellants  have  been  conducted,  and  the  results  are 
compared  with  the  data  obtained  by  Miller  and  Cohen  [13]  for  model  validation. 

METHOD  OF  APPROACH 

DESCRIPTION  OF  PHYSICAL  EVENTS 


In  an  HFCI  experiment,  a  hot  metal  particle  comes  Into  contact  with  a  cold  propellant  at  Time  0. 
In  the  early  phase  of  the  process  (Time  Period  I),  heat  Is  conducted  from  the  hot  particle  to  the 
propellant  without  any  phase  change  or  pyrolysis  (see  Fig.  1).  As  time  progresses,  the  temperature 
of  the  propellant  Increases  and  that  of  the  particle  decreases.  Following  a  period  of  Inert 
heating,  the  propellant  starts  to  decompose,  melt,  and/or  gasify  (Time  Period  II).  During  Time 
Period  II,  a  foam  or  melt  layer  usually  exists  near  the  Interface  region  between  the  hot  fragment 
and  solid  propellant.  Since  the  density  of  the  metal  particle  Is  much  higher  than  the  density  of 
the  decomposed  propellant,  It  displaces  the  decomposed  propellant  and  becomes  imbedded  In  (sinks 
into)  the  propellant,  as  shown  In  Fig.  1.  The  decomposed  species  can  further  react  exothermically 
in  the  gas  phase,  foam  layer,  and/or  condensed  (solid,  liquid)  phase  to  cause  ignition.  If  the  rate 
of  heat  release  in  the  foam  layer  Is  lower  than  that  of  the  heat  loss  to  the  surrounding  and  ambient 
materials,  the  spall  particle  will  be  quenched  without  Introducing  Ignition.  Self-sus;.aihed 
ignition  will  occur  only  If  the  rate  of  heat  generated  by  the  exothermic  reactions  exceeor  the  rate 
of  heat  loss.  In  view  of  the  Importance  of  chemical  reactions  in  the  foam  layer  and  surface  region 
of  solid  propellants,  effective  kinetic  mechanisms  must  be  determined  and  accurate  kinetic  data 
obt  lined.  This  information  and  data  are  then  fed  into  the  theoretical  model  for  realistic 
stmulatlons. 
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TIME  KMtOO  I  (lerow  TMB  nMUTKM  OF  *  MELT  UTBII 


TIME  FCMOO  n  (AFTER  THE  FORMATION  OF  A  MELT  LAYER) 


Fig.  1  Phyaieal  Cvont  of  Hot  FrRgamt 

Conductlv*  Ignition  (HFCl)  ProooaoM 
Roprosantod  by  Various  Casas  In  Two 
niffarant  tins  Parlods 


DETERMINATION  OF  EFFECTIVE  IGNITION  KINETICS  FOR  THE  SAMPLE  PROPELLAMT 

Efforts  to  Identify  alenentary  raactlona  and  neasura  thalr  raaetlon  rates  ware  ruled  out  ai'< 
probably  Infeasible  from  a  technical  standpoint. and  certainly  Inappropriate  to  the  scope  of  the 
modeling  effort  and  the  resources  available  for  nuaerleal  computations. 

Since  the  main  Interest  In  these  reactions  Is  the  energy  source  (or  sink)  terms  in  heat 
transfer  equations,  Differential  Scanning  Calorimetry  (DSC)  was  ohosen  as  a  suitable  teohnlque  for 
measuring  reactive  heat  exchange  of  decomposing  propellant  In  oontaot  with  a  metal  surface. 

Operated  In  ramp  mode,  this  Instrument  can  Increase  the  sample  temperature  linearly  with  time  at 
rates  up  to  100*C/mln.  Although  Ignition  by  hot  fragments  could  Involve  much  higher  Instantaneous 
heating  rates,  these  high  rates  cannot  be  sustained  for  any  appreciable  time  due  to  conductive  and 
convective  loss  mechanisms.  Because  observations  of  hot  fragment  Ignition  under  controlled 
conditions  Indicate  that  5  seconds  or  more  are  required  to  establish  Ignition,  the  lOO*C/mlnute 
limitation  may  not  be  too  restrictive. 

The  propellant  used  In  this  study  has  a  unlnodal  distribution  of  RDX  particle  sixes  (about  5 
micron  average).  The  DSC  teat  samples  were  mlorotomed  to  a  uniform  thickness  of  about  0.4  mm  with 
mass  of  about  1  mg  (t5%)  and  placed  In  covered  and  crimped  pans  perforated  In  four  places  with  a 
straight  pin.  The  pan  perforations  allow  for  pressure  release,  while  retaining  the  bulk  of  any 
decomposition  heat  resulting  from  reactions  occurring  at  or  very  near  the  propellant  surface.  Thus, 
the  technique  does  not  distinguish  between  energetic  reactions  In  the  solid,  liquid,  or  gas  phase; 
the  goal  Is  only  to  measure  net  "localized"  energy  release  (or  absorption). 

DSC  thermograms  for  the  test  propellant  used  here  typically  exhibited  an  endotherm  of  about  20 
cal/g,  starting  at  about  185*C  and  followed  by  an  exotherm  of  about  300  cal/g  wMc  .  peaks  at  about 
260*C.  At  the  onset  of  this  study,  the  Intention  was  to  treat  the  exotherm  as  a  simple  global 
reaction  and  the  endothem  as  a  phase  change.  Ultimately,  both  were  described  as  separate  single 
reactions,  each  with  Its  own  set  of  kinetic  parameters.  A  detailed  analysis  was  reported  In  Ref.  8. 
The  reaction  mechanism  proposed  by  Miller  et  al.  [8]  Is  a  two-step  sequence  of  global  reactions 
represented  by  the  following  equations. 

Solid  ^  Liquid  **2  Gaseous 

Propellant  Melt  *  Products 

where  the  specific  reaction  rate  constants  k^  and  k2  can  be  expressed  In  the  Arrhenius  form 

E, 

-  A^  exp(-  ■pj:)  (2) 


“2  •  *2  rT> 


(3) 
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Ii  the  curve-fitting  process.  Miller  et  al.  found  that  the  Arrhenius  coefficient,  can  be 
satisfactorily  represented  by  a  constant,  and  Aj  can  either  be  approximated  by  constants  or 
expressed  as  a  power  law  of  heating  rate  (r)  in  their  DSC  experiments,  l.e., 

A2  -  constant  or  -  ajf  (**) 


A  typical  comparison  of  the  data  with  the  fitted  curve  is  shown  in  Fig.  2.  The  constant  A2  fitted 
curve  is  slightly  different  from  the  curve  shown  in  Fig.  2,  and  is  equally  acceptable  for  numerical 
simulation  of  HFCI  processes.  Considering  the  enormous  simplification  of  the  actual  chemistry 
afforded  by  the  Irreversible  two-step  idealization,  the  quality  of  representation  is  quits  good. 
Good  representations  of  the  data  are  found  for  a  wide  range  of  heating  rates  and  inert  gas  purge 
flow  rates.  The  DSC  data  for  the  sample  propellant  is  summarized  in  Table  1 ,  where  and  Q;  are 
tne  specific  heat  release  for  endothermic  and  exothermic  reactions. 
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Fig.  2  Comparison  of  Composite  DSC  Thermograms 
(Points)  to  Fitted  Curve  Using  Heating 
Rate  Dependent  Aj  Faotor 


Table  1. 

Kinetic 

Constants 

Determined  from  DSC 

*1 

*2 

El 

(Kcal/ 

E2 

(Kcal/ 

“2 

(s*’*"Vk‘’*) 

»2 

Qi 

QZ 

(s-’) 

(s-») 

mol) 

mol) 

(“) 

(cal/g) 

(eal/j 

1.31x10^’ 

t.98xI0]J[ 

±0.80x10 

69,1 

38.2 

.... 

-21.1 

297 

iO.STxlO-” 

±0.4 

±0.4 

7.79x10!^ 

±2.89x10  ' 

.... 

±5.6 

±26 

*«««  ' 

47.2 

-0.516 

297 

— — 

— — 

— — 

±0.4 

±0.014 

.... 

±26 

Inejrporation  of  Two-Step  Global  Reaction  Kinetics  into  the  HFCI  Theoretical  Model 

Before  describing  the  specific  steps  taken  in  incorporating  the  two-step  global  reaction 
kinetics  into  the  HFCI  model,  it  la  useful  to  give  a  brief  summary  of  the  model  structure.  As 
described  in  Ref.  12,  the  physical  model  is  divided  into  various  regions  (see  Fig.  3):  D  spall 
particle  region;  2)  foam  layer  region  1;  3)  foam  layer  region  2;  and  A)  LOI/A  propellant  region. 

Saon  region  has  its  own  governing  equations  and  boundary  conditions. 

In  order  to  Incorporate  the  two-step  global  kinetics  into  the  HFCI  model,  the  main  modification 
-lonoenirated  on  the  formulation  of  foam  layer  Region  1,  which  is  the  area  controlling  the  onset 
of  Ignition. 

Due  to  the  limited  kinetic  data  resources  available  for  LOVA  propellants  at  the  time  the  HFCI 
model  was  initially  formulated,  a  more  general  consecutive  three-step  reaction  from  solid  propellant 
to  gaseous  products  was  adopted,  l.e., 

**1  ^2 

Solid  — »  Liquid  — »  Gaseous  — ^  Gaseous  (5) 

Propellant  Melt  Reactants  Products 
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Witn  the  conelderatlon  of  the  third  reaction,  the  heat  released  In  the  gas-phase  reaction  could 
lnc''eaae  the  temperature  of  the  gases  in  the  foam  layer.  However,  from  the  DSC  experiments  [S],  a 
two-step  global  reaction  was  proposed,  as  discussed  above. 

Based  upon  two-step  reaction  kinetics  and  energy  flux  balance  shown  in  Fig.  tia,  the  equation 
describing  the  rate  of  change  of  gas-phase  temperature  can  be  derived  as  follows. 


2  “’’a  •  2  • 

'*^s^m^g^pg  Dt  "  ^gas'*’*’*a  ”  ^gas-llq*l^ 


6H°  ,)  -  (C  T*  ♦  nAH®  .)] 


*  "ref''g''lrl2’^s‘'m*^ 


(6) 


wh-ere  n  can  be  considered  as  the  combustion  efficiency  in  the  foam  layer. 


Pig.  iia  Energy  Fluxes  Associated  with  the  Gas 

Phase  in  Region  1  (The  Thickness  of  Foam 
Layer  Is  Highly  Exaggerated) 


SPAU.  MSTICLE 


a. 


Fig.  tb  Energy  Fluxes  Associated  with  the  Liquid 
Phase  In  Region  I  (Thickness  of  Foam 
Layer  Is  Highly  Exaggerated) 


For  the  liquid  phase,  the  balance  of  energy  fluxes  are  given  in  Fig.  Ub,  and  the  equation  of 
liquid-phase  temperature  can  be  written  as 


2  2 

(1-»)vr,L„PjC^j  OT  ■  Vpr«'^pt<^  '  V’ 


♦  (t-<-)vr,  - 


9  gas* 


gas-liq  i 


(7) 


From  the  mass  balances  for  both  gas  and  liquid  phases,  the  following  equations  can  be  obtained  after 
certain  manipulations. 


Dj, 


(1-»)R 


I  •  2  * 

|q"  eitr  -  q" 


(l-rf*)vr2L^Cp^P  ” 


net. 


(1-4')p,*f<»e  “net, 

^(1-<i*fi|i)p 


f<ip 


£C- 


g 


L  vr%/ 
n  SI 


(8) 


where 


«L  •  ■  ‘Sg"g  * 


(9) 
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I  B  r*-  2  • 

V  -1 -  Cq"  i(nrr  -  q«  A,  Q,  ] 

irt  gaa  a  gaa-llq  1,  8„,  L 


Pj-P^ 


8„ 


'net .  r  p 

Fo —  IwiTr - 2^  ''b^^  C(1-i|r»f<i) 


(10) 


In  addition  to  the  above  four  governing  equations  (three  ODES,  one  algebraic  equation),  the 
perfect  gas  law  la  used  to  evaluate  gas  density.  The  energy  flux  balamce  at  the  Interface  between 
liquid  Belt  and  propellant  surface  can  be  given  as 


o’*  •  fl**  *  Q  .  (v  )d 

^llq-solld  ^pr  melt  a  ‘^pr 


(11) 


where  the  heat  of  reaction  Qmelt  defined  as 


Belt 


H® 

“f.l 


♦  C_T 


pt  melt 


)  -  (AH 


T,pr 


(12) 


The  endothermic  reaction  froa  solid  propellant  to  liquid  melt  absorbs  heat  from  liquid  phase 
via  the  last  term  In  Eq.  (11).  The  liquid  phase  Is  heated  by  energy  transferred  from  gas  phase 
through  bubble  Interfaces  In  the  foam  layer.  This  mechanism  Is  modeled  as  dgas-Uq'^i  In  Eq. 

(7).  The  exothermic  reaction  (from  liquid  melt  to  gaseous  products)  releases  heat  to  elevate 
gasophase  temperature  via  the  source  term  associated  with  mg^,^  In  Eq.  (6). 


In  order  to  close  the  above  system,  the  parameter  f,  which  relates  v,  .  (the  radial  velocity 
of  the  liquid  phase  at  the  lateral  surface  of  foam  layer  In  region  1)  and^v^pg  (v^pg  •  f  was 

expressed  as 


f 


C  C- 


«net. 


vr  (-V  ♦€)p_ 
a  m  g 


(13) 


The  above  equation  Is  based  upon  various  limiting  conditions  given  below. 


Case  1 :  f  *  0,  when  the  volumetric  gas  generation 
rate  is  much  smaller  than  the  volumetric 
regression  rate  of  the  melt  layer,  l.e. 


*  ■  «net 


Ca5«  2:  f  ♦  1,  when  wr  (-v  )  £  la  /p 
*  "  *net,  « 


Case  3:  f  »  1,  when  nr  ^(-v  )  <<  m  /p 
*  "  ®net,  .« 

To  avoid  singularity  at  the  onset  of  melting  or  liquefaction  of  LOVA  propellants,  a  small  parameter, 
c.  In  the  order  of  0.0001  is  inserted  In  the  denominator  of  Eq.  (13).  The  coefficient  C  and 
exponent  n  are  taken  to  be  the  unit  in  the  HFCI  simulation. 


The  foam  layer  region  1  thus  contains  six  major  unknown  parameters  (i|t,  Tg,  v  ,f 
Pg,  solved  from  Eqs.  (6),  (7),  (8),  (10),  and  (13).  and  the  perfect  gas  law.  '’*• 


EMPIRICAL  CORRELATIONS  AMD  INPUT  PARAMETERS 


The  exothermic  chemical  reaction  plays  an  Important  role  in  the  HFCI  process,  since  high  rates 
of  heat  release  could  lead  to  runaway  Ignition.  Therefore,  it  is  necessary  to  acquire  an 
appropriate  correlation  for  the  chemical  reaction  rate.  The  kinetic  data  used  in  numerical 
calculations  was  obtained  from  Miller's  DSC  experiments  [8],  The  rate  of  production  of  gaseous  nus- 
In  Region  i  of  tne  foam  layer  is  expressed  as 
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where  A],  Is  the  sane  as  A2  given  in  Table  1.  and  Cg  is  a  constant  sinaller  than  1.  The  temperature 
distribution  In  the  liquid  phase  of  the  foam  layer  Is  not  strictly  uniform.  Most  of  the  liquid  Is 
at  a  temperature  lower  than  the  liquid  surface  temperature  of  bubbles.  This  can  Influence  the  value 
of  Cg.  Also,  the  reacting  surface  layer  of  bubbles  could  be  a  small  fraction  of  the  liquid  in  the 
foam  layer.  Furthermore,  the  Arrhenius  law  Is  highly  nonlinear.  All  of  these  factors  could 
influence  the  effective  mass  generation  rate.  In  the  present  lumped  parameter  analysis  for  the  foam 
layer,  the  selected  value  of  Cg,  In  HPCl  simulation  is  0.02  for  all  cases  studied. 


In  addition  to  the  kinetic  correlation,  suitable  empirical  correlations  for  heat  fluxes  (e.g., 
'5ilq-solld*  ’iliq*  ^gas*  *0^  ^gas-llq^  *1®®  needed  to  properly  predict  heat  transfer  rates  and 
system  eigenvalues,  t^t  Is,  the  sinking  velocity  (v,  or  -Vg)  and^the  propellant  surface  temperature 
(Tpa).  The  heat  transfer  coefficients  for  qii<)-aolld'  Qliq'  Qgas  adopted  In  the  combined 
form  of  forced  convection  and  conduction.  In  the  early  stage  of  time  period  II  (described  In 
physical  events)  or  before  propellant  Ignition,  the  radial  velocity  In  Region  1  of  the  foam  layer 
can  be  quite  high.  Therefore,  the  forced  convection  effect  dominates  heat  transfer  among  the  spall 
particle,  foam  layer,  and  propellant.  However,  In  the  quenched  cases,  the  conduction  heat  transfer 
Is  not  the  main  mechanism  during  most  of  time  period  II.  Some  of  the  heat  transfer  correlations 
used  are  given  In  Ref.  12.  In  consideration  of  various  forma  of  heat  transfer  correlation.  It  is 
noted  that  the  heat  transfer  measurements  available  for  foam  layers  are  very  scarce.  In  fact,  heat 
transfer  rates  of  foams  generated  by  liquefaction  and  pyrolysis  of  solid  propellant  are  completely 
unavailable.  This  la  an  area  requiring  further  research.  The  heat  fluxes  mentioned  above  can  be 
expressed  as  follows: 
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The  heat  flux  at  propellant  surface  (qpr)  Is  calculated  from  the  following  equation,  as 
described  In  Ref.  12. 
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where  the  coefficient  3Ap  results  from  a  third-order  polynomial  approximation  for  the  temperature 
profile.  The  time  variation  of  the  thermal  wave  penetration  depth,  ipp.  Is  determined  from 
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After  5pp  Is  solved  from  Eq.  (20),  the  propellant  surface  temperature  can  be  determined  from  Eq. 
(11)  through  the  use  of  Eqs,  (17)  and  (19). 


The  sinking  velocity  can  then  be  calculated  by 


V  .  (!iiq-S0U8..'..V)  (21) 

* 

In  addition  to  the  above  correlations  and  data  provided  in  Table  1 ,  a  set  of  Input  data  was 
prepared  and  listed  In  Table  2.  The  thermodynamic  and  transport  properties  were  obtained  from  open 
literature  [1ii-19]. 


Ol/ERALL  STRUCTURE  OF  THEORETICAL  MODEL 


The  theoretical  model  consists  of  governing  equations  and  their  associated  Initial  and  boundary 
I'^ndltlons  for  the  hot  particle,  the  propellant,  and  the  foam  layer.  Governing  equations  for  the 
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Table  2. 

Input  Data  Used  In  the  HFCI  Simulation 

Raglon 

Paraaatar 

Valua 

Unit 

Spall 

Pa 

0.3175  X  10"2  or 

D 

0;M851  X  10“2 

m 

Partlcla 

0.3937  X  10"2  or 

m 

0;98  X  10"2 

m 

16.246 

W/m-K 

Cs 

5.02  X  102 

J/kg-K 

Ps 

8;0  X  103 

kg/a3 

2.094  X  10"! 

W/m-k 

1;60  X  103 

kg/n3 

cr 

1;465  X  103 

J/kg-K 

Foaa 

H 

0;832  X  10“3 

kg/m-s 

Layar 

®a  Jg 

l;597  X  io5 

J/oola 

*2. 

l;98  X  10'“ 

1/s 

AKf.i 

*6.695  X  io5 

J/kg 

*^  f  a 

-1.912  X  10“ 

J/kg 

36; 29  X  10“3 

H/m-K 

Cpg 

1.465  X  103 

J/kg-K 

n 

0;80 

^pr 

3.852  X  102  ♦  2.598  »  T(K) 

J/kg-k 

Ppr 

1;678  X  103 

kg/a3 

LOVA 

IpB 

2;094  X  10”' 

H/o-K 

Propallant 

“^,pr 

-7.579  X  105 

J/kg 

^•,Pl 

2.9  X  105 

J/mola 

*1 

1.31  X  1031 

1/s 
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Instantanaous  valuas  of  porosity,  bulk  taaparaturas  of  gas 
and  liquid  phasas,  radial  valoclty  of  foaa  layar,  danslty  of 
tha  gas  phasa,  and  sinking  valoclty  ara  solvad  froa  a  sat  of  mass 
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Flgura  5.  Block  Diagram  Showing  tha  Hathanatlcal  Fomulatlon 

and  Coupling  Rolatlonship  Batwaan  Hot  Partlcla,  Foaa 
Layar,  and  Solid  Propallant 
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hoi  particle  and  the  propellant  are  tranalent  heat  conduction  equations  written  In  two-dloensional 
cylindrical  form.  These  equations  were  recast  Into  ordinary  differential  equations  by  Integral 
methods.  In  these  two  regions,  instantaneous  teoperature  distributions  au*e  solved  from  the 
governing  equations  coupled  to  the  foam  layer  through  flux  balances  at  their  boundaries  (see  Fig.  5). 
The  solution  of  major  unknowns  In  the  foam  layer  la  also  delineated  In  the  same  figure. 


DISCUSSION  OF  RESULTS 

In  the  following,  one  set  of  HFCI  simulation  results  Is  compared  with  experim«ital  results,  and 
two  sets  of  parametric  runs  are  presented  and  discussed  in  detail.  To  study  the  Initial  temperature 
effect  of  the  spall  particle,  two  different  temperatures  (768  and  1000  K)  were  considered  for  the 
same  spall  particle  size  (L5  •  0.3937  cm,  Rg  •  O.3175  cm).  The  calculated  results  from  the  case 
with  an  Initial  temperature  of  768  K  are  compared  with  experimental  data  obtained  by  Miller  and 
Cohen  [13].  The  third  set  of  results  was  obtained  for  studying  the  mass  effect  on  the  HFCI  process 
using  a  large  spall  particle  (Lg  •  0.98  cm,  rg  •  0.1)851  cm). 

Figure  6  shows  a  set  of  propellant  subsurface  temperature  profiles  at  different  times  before 
the  spall  particle  begins  to  sink  Into  the  LOVA  propellant.  As  one  can  see  from  the  temperature 
profile  variation,  the  thermal  wave  penetration  depth  becomes  deeper  as  time  increases.  However, 
the  surface  temperature  remains  approximately  at  the  735  K  range. 

The  axial  temperature  variation  at  the  center  line  of  the  spall  particle  is  shown  in  Fig.  7  for 
different  times  In  time  period  I.  This  figure  Indicates  the  decay  of  maximum  temperature  at  a 
distance  far  from  the  spall  particle  base.  The  temperature  gradient  at  the  base  of  the  spall  drops 
from  a  large  value  upon  Initial  contact,  to  a  much  lower  level  at  1.3  s.  At  the  end  of  time  period 
I,  the  axial  temperature  distribution  becomes  quite  uniform. 


Pig.  6  Calculated  Temperature  Distributions  in  Fig.  7  Calculated  Axial  Temperature 

the  LOVA  Propellant  at  Various  Times  Distributions  Along  the  Centerline  of 

the  Spall  Particle  at  Various  Times 
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In  tine  period  II,  the  spall  particle  sinks  into  the  LOVA  propellant.  The  calculated  sinking 
distance  versus  time  Is  compared  with  measured  data  In  Pig.  8.  The  calculated  sinking  velocity  Is 
also  shown  In  this  figure.  Note  that  the  sinking  velocity  (the  slope  of  the  curve)  increases 
rapidly  at  the  very  beginning,  starts  to  decrease  at  t  <•  0.75  s,  and  approaches  zero  at  the  end  of 
the  calculation  (t-B**  a).  The  final  calculated  sinking  distance  la  about  155  larger  than  the 
measured  one.  This  discrepancy  could  be  caused  by  errors  In  experimental  measurement  as  well  as  any 
Inappropriate  of  heat  transfer  correlation  used  In  the  foam  layer  of  the  theoretical  model.  In 
general,  this  comparison  shows  a  reasonable  agreement  between  calculated  and  measured  results. 

The  calculated  time  variation  of  the  void  fraction  in  the  foam  layer  [(i(t)]  Is  shown  In  Fig.  9. 
When  the  spall  particle  starts  to  sink  Into  the  LOVA  propellant,  Increases  drastically  due  to 
Intensive  heating,  and  then  slows  down  significantly  as  time  Increases  at  the  end  of  computation; 
only  72t  of  the  space  In  the  foam  layer  la  occupied  by  the  gas.  The  calculated  time  variations  of 
center  and  base  temperatures  of  the  spall  particle  are  compared  with  thermocouple  measurements  at 
corresponding  locations  In  Figs.  10  and  11,  respectively.  Comparison  of  the  predicted  base 
temperature  Is  better  than  comparison  of  the  predicted  center  temperature  with  experimental  data. 
Both  of  these  temperatures  decline  monotonlcally  as  the  particle  begins  to  sink  Into  the  LOVA 
propellant  and  approach  asymptotic  levels  near  the  end  of  computation. 

The  calculated  time  variation  of  the  Instantaneous  subsurface  temperature  of  the  first 
thermocouple  location  (6  mm  below  the  Initial  propellant  surface)  compares  well  with  measured  data, 
as  shown  In  Fig.  12.  Like  the  measured  data,  the  calculated  subsurface  temperature  at  the  second 
thermocouple  location  (12  mm  below  the  surface)  shows  no  reponse  at  the  end  of  calculation  (t-3<i  s). 

Comparison  of  bulk  temperatures  of  gas  and  liquid  phases  In  the  foam  layer  la  shown  In  Fig.  13. 
During  time  period  II,  the  gas-phase  temperature  Is  always  higher  than  that  of  the  liquid.  The 
temperature  difference  decreases  monotonlcally  and  reaches  a  small  difference  near  the  end  of 
computation. 


TiMe,  1-1,.  » 

Fig.  8  Comparison  of  Calculated  Time  Variations 
of  Sinking  Distances  with  Experimental 
Data  and  Calculated  Time  Variation  of 
Sinking  Velocity  of  the  Spall  Particle 
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Fig.  9  Calculated  Time  Variation  of  the  Void 
Fraction  In  the  Foam  Layer 
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Fig.  10  Conpsrlson  of  Caloulatsd  Taaparatura- 
Tl«a  Traea  at  tiia  Cantor  of  tha  Spall 
Parttela  with  Exporlaontal  Data 
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Fig.  11  Coaparlaon  of  tha  Taaperatura-Tlae  Traea 
at  tha  Bottoa  Surfaca  of  tha  Spall 
Partlola  with  Exparlnantal  Data 


Fig.  12  Conparlaona  of  Prodlotod  Taaparatura- 
Tlna  Traeaa  at  tha  First  and  Saeond 
Thamocoupla  Locations  (6  and  12  m 
Balow  tha  Propollant  Surfaoo)  with 
Exparlaantal  Data 


Fig.  13  Calculated  Tanperatura-TlRa  Traeaa  of 

Caa  and  Liquid  Phases  In  the  Foa«  Layer 
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In  addition  to  the  quenched  case  discussed  above,  two  separate  HFCI  tests  were  conducted  by 
Miller  and  Cohen  [13],  using  different  spall  particle  sizes  and  Initial  temperature.  The  small 
cylindrical  particle  has  a  radius  of  0.3175  cm  and  a  length  of  .3937  cm,  with  an  Initial  temperature 
of  986  K  at  the  time  of  contact  with  the  LOVA  propellant.  The  large  cylindrical  particle  has  a 
radius  of  0.4651  cm  and  a  length  of  0.98  cm,  with  an  initial  temperature  of  951  K  at  the  time  of 
contact  with  the  LOVA  propellant.  In  both  cases,  the  propellant  ignited.  The  pyrolyzed  gases  are 
not  luminous  near  the  base  region  of  the  spall  particle.  The  observed  luminous  flame  Is  above  the 
Initial  propellant  surface.  The  flame  widens  abruptly  at  5.2  and  1.2  s,  corespor.dlng  to  the  amal' 
and  large  particles,  respectively.  Widening  of  the  flame  could  be  considered  as  onset  of  sustained 
Ignition  from  the  experimental  point  of  view. 

«<.lldate  tne  tneoretioal  .tiouel  witn  the  Ignited  cases,  both  tests  were  simjlsted.  =esults 

are  discussed  below.  Figure  14  shows  the  comparison  of  calculated  trajectories  of  spall  particles 
with  experimental  data.  The  small  particle  has  a  slightly  higher  Initial  temperature  than  the  large 
one,  and  hence  sinks  Into  the  LOVA  propellant  at  a  slightly  higher  velocity  (see  Fig.  15).  At  a 
later  time,  the  velocity  of  the  large  particle  la  higher  since  it  contains  more  thermal  energy.  The 
calculated  trajectory  shown  In  Fig.  14  overpredicts  the  experimental  data  by  25>  at  the  end  of  the 
run.  However,  the  agreement  la  still  acceptable  for  most  of  the  event.  For  these  two  cases,  the 
particles  sink  eontlnoualy  as  the  propellant  reaches  sustained  ignition  conditions.  The  calculated 
sinking  velocity  variations  with  respect  to  time,  shown  In  Fig.  15,  exhibit  spikes  during  the 
Initial  Interval  when  the  particles  are  extremely  hot. 


TIME  ,  l-l,,« 

FIr.  14  Comparison  of  Calculated  Time  Variations 
of  Sinking  Distances  with  Experimental 
Data 
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Fig.  15  Calculated  Time  Variations  of  Sinking 
Velocity 


Figure  16  shows  the  comparison  of  calculated  teraperature-tlme  traces  at  the  base  of  the  spall 
particle  with  experimental  data.  During  the  early  phase,  the  base  temperature  for  small  particles 
Is  higher  than  that  for  large  particles;  the  trend  reverses  at  the  later  phase  of  Time  Period  II 
(when  spall  particles  sink  Into  the  LOVA  propellant).  Calculated  results  are  In  good  agreement  with 
experimental  data  In  terms  of  the  trend  mentioned  above,  as  well  as  the  magnitude  and  slope  of  these 
traces.  It  Is  Interesting  to  note  that  the  average  temperature  at  the  spall  particle  base  drops 
continuously  even  after  onset  of  sustained  ignition;  this  Is  believed  to  be  the  result  of  Incomplete 
combustion  In  the  foam  layer.  Gaseous  pyrolysis  products  In  the  bubbles  take  time  to  reach  the 
fully  reacted  state.  Heat  release  In  the  foam  zone  Is  sufficiently  low  to  produce  a  relatively  c-ol 
layer  of  foam  material  which  reduces  the  energy  of  the  spall  particle.  The  above  effect  can  be  seen 
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from  the  plot  of  the  calculated  bulk  temperaturea  of  the  gaa  and  liquid  phaaes  In  the  foam  layer 
(see  Fig.  17).  The  temperaturea  of  both  phaaea  are  lower  than  the  average  temperature  of  the  apall 
particle  base .  The  calculated  void  fraction  variations  with  reapect  to  time  la  shown  in  Fig.  18. 
After  the  Initial  Intensive  heating  In  time  period  II,  the  void  fraction  reduces  slightly  with 
respect  to  tisw. 


Comparisons  of  calculated  temperature-time  traces  at  the  center  of  spall  particles  with 
erperimental  data  la  shown  In  Fig.  19.  Agreement  la  quite  reasonable,  and  the  trend  is  Identical  to 
that  described  for  spall  particle  base.  Before  completion  of  the  computation  of  10  a  for  the  time 
oeriod  II,  the  center  temperaturea  are  always  :'.<gner  than  those  at  the  spall  particle  base. 
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Fig.  16  Caap«rl8on  of  Calculated  Tamparatura- 
Tiaa  Traoaa  at  tha  Baaa  of  tha  Spall 
Fartlola  with  Cxparlsantal  Data 
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"ig.  17  Calculated  Bulk  Tamparaturas  of  tha  Gaa 
and  Liquid  Phases  in  tha  Foan  Layer 
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Fig.  18  CmleulatAd  Tin*  VmrlAtlona  of  Void 
Fractions  In  Foaa  Layars 
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Fig.  19  CoRparlson  of  Calculated  Teaperature- 
Tlae  Traces  at  the  Center  of  Spal. 
Fartlcles  with  Cxperlnental  Data 


Based  upon  these  comparisons.  It  Is  quite  obvious  that  the  HFCI  model  Is  able  to  predict  the 
quenched  and  Ignited  cases  using  different  spall  particle  sizes.  Also,  by  comparing  results  of  the 
small  particle  (L  •  0.3937  cm,  r,  -  0.3175  cm)  at  an  Initial  temperature  of  986  K  with  those  for  the 
same  sized  particle  at  1000  K,  one  can  observe  the  strong  difference  In  solution  from  a  quenched 
case  to  a  runaway  Ignition.  According  to  the  test  data  of  Miller  and  Cohen  [13],  the  small  particle 
can  Ignite  the  LOVA  propellant  at  986  K,  which  Is  m  agreement  with  the  prediction  of  ignition  at 
986  K. 

SUMHART  AKD  CONCLUSIONS 

1)  The  two-step  chemical  kinetic  mechanism,  proposed  by  Hiller  et  al.  [8],  based  upon  DSC 
measurements,  was  Incorporated  Into  the  theoretical  model  [12]  for  simulating  hot  fragment 
conductive  Ignition  processes  of  LOVA  propellants. 

2)  Calculated  results  compared  well  with  experimental  data  In  temperature-time  traces  of  spall 
particle  and  propellant. 

3)  Reasonable  agreement  was  achieved  between  theoretical  prediction  and  experimental  measurements 
In  terms  of  the  terminal  position  and  particle  trajectory  of  the  quenched  spall  particle  which 
Is  partially  submerged  in  the  LOVA  propellant.  The  predicted  trajectory  of  the  spall  particle 
also  agrees  quite  well  with  measured  data  for  Ignited  cases. 

4)  The  controlling  mechanism  In  the  HFCI  process  Is  the  competition  between  the  rate  of  heat 
release  from  exothermic  reactions  Introduced  by  the  spall  particle  and  the  rate  of  heat  loss  to 
the  surrounding. 

5)  The  present  two-step  reaction  model  gave  areasonable  numerical  solution  In  HFCI  simulation. 
However,  a  more  complicated  reaction  model  may  be  needed  for  various  types  of  LOVA  propellants. 
Also,  additional  heat  transfer  correlations  in  the  foam  zone  of  the  pyrolysis  product  of  various 
propellants  should  be  developed  for  more  accurate  predictions  of  HFCI  processes. 

6)  From  HFCI  point  of  view,  it  la  highly  desirable  to  have  LOVA  propellants  with  cool  pyrolysis 
products  In  the  foam  layer  which  can  effectively  quench  the  hot  spall  particle.  Also,  It  is 
desirable  to  have  LOVA  propellants  which  require  extended  distances  for  development  of 
flames . 
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Validation  of  a  Theoretical  Model 
for  Hot  Fragment  Conductive  Ignition  Processes 
of  LOVA  Propellants 
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A  comprehensive  theoretical  model  and  an  efficient  numerical 
program  have  been  developed  to  simulate  the  hot  fragment  conductive 
Ignition  (HFCI)  processes  for  characterizing  the  degree  of 
vulnerability  of  various  gun  propellants.  In  the  formulation  of  the 
theoretical  model,  three  seta  of  governing  equations  and  their 
boundary  conditions  were  derived  for  different  regions  consisting  of 
the  hot  spall  particle,  the  LOVA  propellant,  and  the  foam  layer 
produced  by  liquefaction,  pyrolysis,  and  decomposition  of  the 
propellant.  In  terms  of  chemical  kinetic  scheme,  a  two-step  reaction 
model  Is  proposed,  baaed  upon  DSC  experiments.  The  model  contains  an 
endothermic  decomposition  process  followed  by  an  exothermic  pyrolysis 
of  the  propellant.  Calculated  results  compared  well  with  experimental 
data  In  temperature-time  traces,  trajectory  of  spall  rt'  .lules  sinking 
Into  the  LOVA  propellant,  and  terminal  position  of  the  particle. 


NOMENCLATURE 

A  Arrhenius  frequency  factor  of  chemical 

reaction  of  melt  or  foam  layer  in 
Region  1,  kg/n3-s 

A]  Interface  area  between  liquid  melt  and 

gas  bubbles,  m^ 

Cp  constant-pressure  specific  heat, 

J/kg-K 

C,  specific  heat  of  spall  particle, 

J/kg-K 

Jp  averaged  diameter  of  bubbles  generated 
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from  gasification  process,  m 
Ea  activation  energy  of  chemical  reaction 

of  melt  or  foam  layer,  J/mole 
heat  of  formation,  J/kg 

l<m  height  of  melt  or  foam  layer  In  Region 

1 ,  m 


■gnet, 

Ql 


^^raelt 

*^gas 

»  It 

'^gas-llq 


e" 

<111q 

^lllq-solld 


net  rate  of  production  of  gaseous  mass 
In  Region  1  of  foam  layer,  kg/s 
heat  release  due  to  exothermic 
reaction  from  liquid  melt  to  gaseous 
products,  J/kg 

heat  of  endothermic  reaction  from 
solid  propellant  to  liquid  melt,  J/kg 
heat  flux  transferred  to  gas  phase 
from  spall  particle  base,  W/m^ 
heat  flux  transferred  to  liquid 
phase  from  gas  phase  through 
Interfaces,  W/m^ 
heat  flux  transferred  to  liquid 
phase  from  spall  particle-  base,  W/ra^ 
heat  flux  transferred  to  propellant 
surface  from  foam  layer,  W/ra2 
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Qpp 


<lspall 

Ru 

'*3 

Tb(Tsat) 

Is 

^melt 

Tp 

Tp3 

Ts 

''ir 


Vs  O'"  -Vm 


z 

z* 


heat  flux  tpansferred  to  LOVA 
propellant  from  melt  or  foam  layer  In 
Region  1,  W/m^ 

heat  flux  at  base  of  spall  particle, 
W/m2 

Universal  gas  constant, 

J/kg-mole  K 

radius  of  spall  particle,  m 
boiling  temperature  of  liquid  melt,  K 
gas  temperature,  K 
bulk  temperature  of  melt  or  foam 
layer,  K 

melting  temperature  of  LOVA 
propellant,  K 

temperature  of  LOVA  propellant,  K 
average  surface  temperature  of  LOVA 
propellant  exposed  to  air,  K 
temperature  of  spall  particle,  K 
radial  velocity  of  melt  or  foam 
layer  In  Region  1,  m/s 
sinking  velocity  of  spall  particle, 
rn/3 

axial  distance  above  base  of  spall 
particle,  a 

axial  distance  below  Instantaneous 
surface  of  propellant  under  foam 
layer,  m 


Greek  Symbols 

a  thermal  dlffuslvlty,  m^/s 

4pp  thermal  penetration  depth  from 

propellant  surface,  m 
l  thermal  conductivity,  W/m-K 

p  density,  kg/m3 

y  average  porosity  (void  fraction)  of 

foam  layer  In  Region  1 
u  dynamic  viscosity,  kg/m-s 


Subscripts 


g 

I 

P 

melt 
pr  or  p 


gas-phase 
Initial  condition 
liquid-phase 
melt  liquid  melt 
propellant 


In  the  selection  of  low  vulnerability 
ammunition  (LOVA)  propellants,  go/no-go  Ignition 
tests  with  hot  metallic  elements  are  often 
conducted  to  simulate  hot  fragment  conductive 
Ignition  (HFCI)  processes.  These  hot  fragments 
can  be  considered  as  those  generated  from  the 
penetration  of  armor  plates  by  shaped  charge  Jets 
or  kinetic  energy  penetrators.  One  way  to  negate 
the  threat  of  propellant  Ignition  by  hot  spalls  Is 
to  use  a  propellant  which  Is  resistant  to 
conductive  Ignition.  A  number  of  investigations 
on  this  subject  have  been  conducted,  both 
theoretically  and  experimentally,  in  recent  years 
[1-12].  Since  a  review  of  pertinent  work  la  given 
In  the  paper  by  Kuo  et  al.  [12],  detailed  reviews 
are  omitted  here.  Only  major  findings  of  HFCI 
processes  are' listed  below. 

1.  Binders  of  nltramlne  composite 
propellants  were  found  to  have  a  strong 
effect  on  conductive  Ignition.  Some 
binders  can  act  as  fire  retardant 
coolants. 

2.  The  relative  susceptibility  of  LOVA 
propellants  to  Ignition  by  spall 
particles  can  be  determiner  from 
Ignition  map  based  upon  the  plot  of  the 
initial  temperature  versus  weight  of  the 
spall  fragment. 

3.  The  contact  resistant  between  hot 
particle  and  propellant  significantly 
affects  the  delay  time  for  onset  of 
runaway  Ignition. 

M.  Binders  which  endothermically  decompose 
under  acid  catalysis  are  more  desirable 
for  LOVA  propellant  binder  ingredients. 

5.  The  Ignition  temperatures  of  LOVA 
propellants  are  similar  to  their 
nltramlne  fillers. 

6.  For  certain  LOVA  propellants, 
liquefaction  and  bubble  formation  were 
observed  during  conductive  heating. 

This  suggests  that  the  thermal 
Insulation  properties  of  binder  products 
may  be  Important  In  the  determination  of 
Ignition  sensitivities  of  LOVA 
propellants. 

7.  Based  upon  surface  thermocouple 
measurements  during  HFCI  tests,  the 
Initial  condensed  phase  reactions  are 
endothecmlc.  The  endothermlclty  of 
these  reactions  showed  good  correlation 
with  propellant  sensitivity  for  thermal 
Ignition. 

8.  The  minimum  spall  particle  temperature 
(for  a  fixed  spall  particle  weight) 
depends  strongly  on  binder  composition. 

9.  Chemical  reactions  leading  to  Ignition 
of  certain  LOVA  propellants  during  HFCI 
experiments  can  be  idealized  as  a 
two-step  sequence  of  global  reactions 
in  which  an  endothermic  reaction  is 
followed  by  an  exothermic  reaction. 


•  room  condition 

1  Region  1  of  foam  layer 

Superscript 


Several  attempts  were  made  to  model  the  hot 
fragment  conductive  Ignition  processes  [9-1']; 
however,  these  models  have  serious  limitations  In 
simulating  actual  HFCI  processes.  The  model 
proposed  in  Ref.  9  la  strictly  one-dlme.islonal , 
and  allows  no  space  for  the  products  In  the  gas  or 
liquid  phases  to  exit  at  the  Interface  between 


T'  •  T  -  Tref 


propellant  aanples  and  hot  Inert  fragments.  The 
model  proposed  In  Ref.  10  requires  no  depletion  of 
propellant  when  the  hot  fragment  penetrates  the 
combustible  material.  Also,  no  surface  reaction 
is  allowed  in  the  model.  An  interesting 
theoretical  analysis  of  the  Ignition  of  reactive 
solids  by  direct  contact  wltl  a  hot  inert  body  was 
performed  by  Linan  and  Klndelan  [11].  Their 
analysis  considered  both  cylindrical  and  spherical 
geometries.  Asymptotic  solutions  for  large 
activation  energies  were  obtained.  The  solution 
Indicated  the  existence  of  an  ignition  boundary 
beyond  which  runaway  Ignition  occurred.  Although 
the  theoretical  results  Indicated  the  correct 
trend  in  terms  of  inert  particle  initial 
temperature  and  particle  size,  their  solutions 
were  not  compared  with  any  experimental  data. 
Furthermore,  the  assumption  of  high  activation 
energies  are  Invalid  for  various  LOVA  propellants. 
There  was  no  provision  for  the  development  of  a 
foam  layer  at  the  Interface  region  between  hot 
fragment  and  propellant. 

A  comprehensive  theoretical  model  was 
recently  developed  by- Kuo  et  al.  [12]  to  describe 
the  heat  transfer  processes  in  hot  fragment  and 
solid  propellant,  as  well  as  the  development  of  a 
foam  and/or  melt  layer  which  has  been 
experimentally  observed  at  the  interface  zone 
between  solid  propellant  surface  and  hot  spall 
particle.  In  this  model,  the  liquefaction, 
surface  reaction,  and  gas-phase  reactions  in  the 
foam  layer  were  considered.  The  effect  of  binder 
composition  can  also  be  studied  by  using  measured 
kinetic  data  from  separate  DSC  experiments  for 
U)VA  propellants  with  various  binders.  Thus,  the 
ability  of  resisting  conductive  ignition  of  new 
propellant  formulations  can  be  evaluated  for 
achieving  reduced  system  vulnerability. 

Although  the  model  presented  in  Ref.  12 
includes  both  reacting  and  non-reacting  cases, 
only  the  non-reacting  part  of  the  model  was 
verified  by  ice-melting  and  evaporation 
experiments.  The  main  objective  of  this  paper  is 
to  verify  the  chemically  reacting  part  of  the 
model.  The  latest  reaction  mechanism  and  data 
obtained  by  Miller  et  al.  [8]  have  been 
Incorporated  into  the  model.  Numerical  simulation 
of  HFCI  processes  of  certain  LOVA  propellants  have 
been  conducted,  and  the  results  are  compared  with 
the  data  obtained  by  Miller  and  Cohen  [131  for 
model  validation. 

METHOD  OF  APPROACH 
Description  of  Physical  Events 

In  an  HFCI  experiment,  a  hot  metal  particle 
comes  into  contact  with  a  cold  propellant  at  Time 
0.  In  the  early  phase  of  the  process  (Time  Period 
I),  heat  is  conducted  from  the  hot  particle  to  the 
propellant  without  any  phase  change  or  pyrolysis 
(see  Fig.  1).  As  time  progresses,  the  temperature 
of  the  propellant  increases  and  that  of  the 
particle  decreases.  Following  a  period  of  inert 
heating,  the  propellant  starts  to  decompose,  melt, 
and/or  gasify  (Time  Period  II).  During  Time 
Period  II,  a  foam  or  melt  layer  usually  exists 
near  the  Interface  region  between  the  hot  fragment 
and  solid  propellant.  Since  the  density  of  the 
metal  particle  is  much  higher  than  the  density  of 
the  decomposed  propellant,  it  displaces  the 
decomposed  propellant  and  becomes  Imbedded  in 
(sinks  into)  the  propellant,  as  shown  in  Fig.  1. 


TIME  PERIOO  I  (  BEFDK  THE  FORMATION  OF  A  MELT  LATER ) 

h 

‘r  I II  '^"**°* 

•  SOLID  PROPELLANT 


TIME  PERIOD  0  (AFTER  THE  FORMATION  OF  A  MELT  LAYER) 

LATER 


Fig.  1  Physical  Event  of  Hot  Fragment 

Conductive  Ignition  (HFCI)  Processes 
Represented  by  Various  Cases  in  Two 
''•fferent  Time  Periods 

The  decomposed  species  can  further  react 
exothermically  in  the  gas  phase,  foam  layer, 
and/or  condensed  (solid,  liquid)  phase  to  cause 
ignition.  If  the  rate  of  heat  release  in  the  foam 
layer  is  lower  than  that  of  the  heat  loss  to  the 
surrounding  and  ambient  materials,  the  spall 
particle  will  be  quenched  without  introducing 
ignition.  Self-sustained  ignition  will  occur  only 
if  the  rate  of  heat  generated  by  the  exothermic 
reactions  exceeds  the  rate  of  heat  loss.  In  view 
of  the  importance  of  chemical  reactions  in  the 
foam  layer  and  surface  region  of  solid 
propellants,  effective  kinetic  mechanisms  must  be 
determined  and  accurate  kinetic  data  obtained. 

This  information  and  data  are  then  fed  into  the 
theoretical  model  for  realistic  simulations. 

Determination  of  Effective  Ignition  Kinetics  for 
the  Sample  Propellant 

Efforts  to  Identify  elementary  reactions  and 
measure  their  reaction  rates  were  ruled  out  as 
probably  Infeasible  from  a  technical  standpoint, 
and  certainly  inappropriate  to  the  scope  of  the 
modeling  effort  and  the  resources  available  for 
numerical  computations. 

Since  the  main  Interest  in  these  reactions  is 
the  energy  source  (or  sink)  terms  in  heat  transfer 
equations.  Differential  Scanning  Calorimetry  (DSC) 
wu  chosen  as  a  suitable  technique  for  measuring 
reactive  heat  exchange  of  decomposing  propellant 
in  contact  with  a  metal  surface.  Operated  in  ramp 
mode,  this  Instrument  can  increase  the  sample 
temperature  linearly  with  time  at  rates  up  to 
100°C/oln.  Although  Ignition  by  hot  fragments 
could  involve  much  higher  instantaneous  heating 
rates,  these  high  rates  cannot  be  sustained  for 
any  appreciable  time  due  to  conductive  and 
convective  loss  mechanisms.  Because  observations 
of  hot  fragment  ignition  under  controlled 
conditions  Indicate  that  5  seconds  or  more  are 
required  to  establish  ignition,  the  1 0CC/ralnute 
limitation  may  not  be  too  restrictive. 

The  propellant  used  in  this  study  has  a 
unimodal  distribution  of  RDX  particle  sizes  (about 
5  micron  average).  The  DSC  test  samples  were 
mlcrotomed  to  a  uniform  thickness  of  about  O.u  mm 
with  mass  of  about  1  mg  (t5S)  and  placed  in 
covered  and  crimped  pans  perforated  in  four  places 
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with  a  straight  pin.  Th«  pan  perforations  allow 
for  pressure  release,  while  retaining  the  bulk  of 
any  decoapositlon  heat  resulting  fron  reactions 
occurring  at  or  very  near  the  propellant  surface. 

Thus,  the  technique  does  not  distinguish  between 
energetic  reactions  in  the  solid,  liquid,  or  gas 
phase;  the  goal  is  only  to  measure  net  "localized" 
energy  release  (or  absorption). 

DSC  thermograms  for  the  teat  propellant  used 
here  typically  exhibited  an  endothern  of  about  20 
cal/g,  starting  at  about  185*C  and  followed  by  an 
exotherm  of  about  300  cal/g  which  peaks  at  about 
260*C.  At  the  onset  of  this  study,  the  intention 
was  to  treat  the  exotherm  as  a  simple  global 
reaction  and  the  endotherm  as  a  phase  change. 
Ultimately,  both  were  described  as  separate  single 
reactions,  each  with  its  own  set  of  kinetic 
parameters.  A  detailed  analysis  was  reported  in 
Ref.  8.  The  reaction  mechanism  proposed  by  Miller 
et  al.  [8]  la  a  two-step  sequence  of  global 
reactions  represented  by  the  following  equations. 

Solid  ‘‘t  ^  Liquid  *^2^  Caseous  .j. 

Propellant  Melt  Products  ' 

where  the  specific  reaction  rata  constants  k<  and 
k2  can  be  expressed  in  the  Arrhenius  form 

E, 

k,  -  A,  exp(-  •^)  (2) 


kj  •  *2 

u 

tn  the  curve-fitting  process,  Hiller  et  al.  found 
that  the  Arrhenius  coefficient.  A;  can  be 
satisfactorily  represented  by  a  constant,  and  Aj 
can  either  be  approximated  by  constants  or 
expressed  as  a  power  law  of  heating  rata  (r)  in 
their  DSC  experiments,  i.e., 

.2 

A^  -  constant  or  A^  •  a^r  (U) 

A  typical  comparison  of  the  data  with  the  fitted 
curve  is  shown  in  Fig.  2,  The  constant  Aj  fitted 
curve  is  slightly  different  from  the  curve  shown 
in  Fig.  2,  and  is  equally  acceptable  for  numerical 
simulation  of  HFCI  processes.  Considering  the 
enormous  simplification  of  the  actual  chemistry 
afforded  by  the  irreversible  two-step 
idealization,  the  quality  of  representation  is 
quite  good.  Good  representations  of  the  data  are 
found  for  a  wide  range  of  heating  rates  and  inert 
gas  purge  flow  rates.  The  DSC  data  for  the  sample 
propellant  is  summarized  in  Table  1. 


where  0;  and  Qj  are  the  specific  heat  release  for 
endothermic  and  exothermic  reactions. 
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Fig.  2  Comparison  of  Composite  DSC  Thermograms 
(Points)  to  Fitted  Curve  Using  Heating 
Rate  Dependent  Aj  Factor 


Incorporation  of  Two-Step  Global  Reaction 
Kinetics  into  the  HFCI  Theoretical  Model 


Before  describing  the  specific  steps  taken  in 
incorporating  the  two-step  global  reaction 
kinetics  into  the  HFcI  model,  it  la  useful  to  give 
a  brief  summary  of  the  model  structure.  As 
described  in  Ref.  12,  the  physical  model  Is 
divided  into  various  regions  (see  Fig.  3):  D 
spall  particle  region;  2)  foam  layer  region  t;  3) 
foam  layer  region  2;  and  A)  LOVA  propellant  region. 
Each  region  has  its  own  governing  equations  and 
boundary  conditions.  •' 

In  order  to  Incorporate  the  two-step  global 
kinetics  into  the  HFCI  model,  the  main 
modification  is  concentrated  on  the  formulation  of 
foam  layer  Region  1,  which  is  the  area  controlling 
the  onset  of  ignition. 

Due  to  the  limited  kinetic  data  resources 
available  for  LOVA  propellants  at  the  time  the 
HFCI  model  was  initially  formulated,  a  more 
general  consecutive  three-step  reaction  from  solid 
propellant  to  gaseous  products  was  adopted,  I.e., 

k  k  k  * 

1  *2  *3 

Solid  — »  Liquid  — ♦  Caseous  Gaseous 

Propellant  Melt  Reactants  Products 

(5) 

With  the  consideration  of  the  third  reaction,  the 
heat  released  in  the  gas-phase  reaction  could 
Increase  the  temperature  of  the  gases  In  the  foam 
layer.  However,  from  the  DSC  experiments  [8],  a 
two-step  global  reaction  was  proposed,  as 
discussed  above. 


Table  1. 

Kinetic 

Constants 

Determined 

from  DSC 

Al 

*2 

El 

Eg 

•2 

bj 

Qi 

02 

(s-’) 

(s-i) 

(Kcal/ 

mol) 

(Kcal/ 

mol) 

(s'’*‘’/k'’*) 

(  — ) 

(cal/g) 

(cal/( 

1.31x10^’ 

±0.57x10-’ 

1.98xlo!* 

±0.80x10 

69.  A 

38.2 

— 

-21.1 

297 

±0.4 

±0.4 

7.79x101, 
±2.89x10  ' 

±5.6 

±26 

—  — 

47.2 

-0.516 

297 

— 

— 

±0.4 

±0.01 4 

±26 
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Baaed  upoa  two-atep  reaction  klnetlca  and 
energy  riux  balance  ahown  In  Fig.  <ia,  the  equation 
deacrlblng  the  rate  of  change  of  gaa-phaae 
teaperature  can  be  derived  aa  folloua. 


2  °^g  •  2  . 

♦’^a^'a.Og'^pg  oT  •  ‘’S.3*’"'a  '  ^Jaa-llq*!, 


*  '’^ref‘'g''lrt2’^3'-ffl^' 


(6) 
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Fig.  4a  Energy  Fluxes  Associated  with  the  Gas 

Phase  In  Region  1  (The  Thickness  of  Foam  • 
Layer  is  Highly  Exaggerated) 


Fig.  4b  E.nergy  Fluxes  Associated  with  the  Liquid 
Phase  in  Region  i  (Thickness  of  Foam 
Layer  Is  Highly  -xaggerated 


where  n  can  be  considered  aa  the  combustion 
efficiency  In  the  foam  layer. 

For  the  liquid  phase,  the  balance  of  energy  fluxes 
are  given  In  Fig.  4b,  and  the  equation  of 
liquid-phase  teaperature  can  be  written  as 


?  2 

*  Qllq  -  qliq,30iid"'s  *  '’Jae-llq*!, 

(7) 

From  the  mass  balances  for  both  gas  and  liquid 
phases,  the  following  equations  can  be  obtained 
after  certain  manipulations. 


P,-P_  net,  r  p 

2^  *  IT  7f  ''a'/ 

tgrasml  pgg 


In  addition  to  the  above  four  governing 
equations  (three  ODES,  one  algebraic  equation), 
the  perfect  gas  law  Is  used  to  evaluate  gas 
density.  The  energy  flux  balance  at  the  Interface 
between  liquid  melt  and  propellant  surface  can  be 
given  as 

*^llq-solld  "  ^pr  *  ®aelt^'^s^*’pr 


0^ 

Dt 


(1-»)R 


(1-*-f*)irrXC„^P 
3  m  pg 


iq"  i»iTr  -  q"  , ,  A, , 
'^gas*^  3  ^gas-llq  11 


RT  ,p  V,  .2xr  L  ijif  ♦  a 
ref*^g  trt  s  g, 


net. 


where  the  heat  of  reaction  Opielt  1®  defined  as 


pr 


Sr  "melt' 


^(l-i)i)pj*fiKPg^ 

M1-V*fiK)Pg 


net. 


f*P 


2L. 


,  (I-4irfi(i)L  p  n 


where 


§0  *0 

a  -  (C  T.  ♦  AH.  -  (C  T  *  nAH,  ) 
L  pi  t  f,l  Pg  g  f.g 


(8) 


(9) 


The  endothermic  reaction  from  solid 
propellant  to  liquid  melt  absorbs  heat  from  liquid 
phase  via  the  last  term  in  Eq.  (ll).  The  liquid 
phase  Is  heated  by  energy  transferred  from  gas 
phase  through  bubble  Interfaces  In  the  foam  layer. 
This  mechanism  Is  modeled  as  9gas-llq’*l  (h  Eq. 
(7).  The  exothermic  reaction  (from  liquid  melt  to 
gaseous  products)  releases  heat  to  elevate 
gas-phase  temperature  via  the  source  term 
associated  with  In  Eq.  (6). 
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In  order  to  close  the  above  systetn,  the 
parameter  t,  which  relates  (the  radial 

velocity  of  the  liquid  phase  at  the  lateral 
surface  of  foam  layer  in  region  1)  and  v,rg  (''erg 

•  f  V.  was  expressed  as 
in 


f 


“net. 


c  C- 


(»3) 


The  above  equation  is  based  upon  various  limiting 
conditions  given  below. 

Case  1 :  f  0,  when  the  volumetric  gas  generation 
rate  la  much  smaller  than  the  volumetric 
regression  rate  of  the  melt  layer,  l.e. 


»r  (-V  )  »  m  /p 

*  “  «net,  « 


Case  2:  f  «  1,  when  «r  {-v  )  a  /p 
*  "  “net,  “ 


Case  3:  f  >>  1 .  when  «  m  /p^ 

“net^ 

To  avoid  singularity  at  the  onset  of  melting  or 
liquefaction  of  LOVA  propellants,  a  small 
parameter,  c.  In  the  order  of  0.0001  is  inserted 
In  the  denominator  of  Eq.  (13)'  The  coefficient  C 
and  exponent  n  are  taken  to  be  the  unit  in  the 
HFCI  simulation. 

The  foam  layer  region  1  thus  contains  six 
major  unknown  parameters  (ip,  Tg,  T^, 

Pg,  solved  from  Eqs.  (6),  (7),  (8),  (107,  and 
(13),  and  the  perfect  gas  law. 

Empirical  Correlations  and  Input  Parameters 

The  exothermic  chemical  reaction  plays  an 
Important  role  In  the  HFCI  process,  since  high 
rates  of  heat  release  could  lead  to  runaway 
Ignition.  Tht.  efore.  It  Is  necessary  to  acquire 
an  appropriate  correlation  for  the  chemical 
reaction  rate.  The  kinetic  data  used  In  numerical 
calculations  was  obtained  from  Miller's  DSC 
experiments  [8].  The  rata  of  production  of 
gaseous  mass  In  Region  1  of  the  foam  layer  Is 
expressed  as 

E 

m  -cl  vr  '(1-<()n,A,exp(-  )  (I*) 

“net,  _  "  ■  »  ‘  “u^l 

where  is  the  same  as  Ag  given  In  Table  i,  and 
c„  Is  a  constant  smaller  than  1 .  The  temperature 
distribution  In  the  liquid  phase  of  the  foam  layer 
Is  not  strictly  uniform.  Host  of  the  liquid  Is  at 
a  temperature  lower  than  the  liquid  surface 
temperature  of  bubbles.  This  can  Influence  the 
value  of  Cj,.  Also,  the  reacting  surface  layer  of 
bubbles  could  be  a  small  fraction  of  the  liquid  In 
the  foam  layer.  Furthermore,  the  Arrhenius  law  Is. 
highly  nonlinear.  All  of  these  factors  could 
Influence  the  effective  mass  generation  rate.  In 
the  present  lumped  parameter  analysis  for  the  foam 
layer,  the  selected  value  of  c„,  in  HFCI  simulation 
is  0.02  for  all  cases  studied. 


In  addition  to  the  kinetic  correlation, 
suitable  empirical  correlations  for, heat  fluxes 

(e.g.,  4Ilq-solld>  "illq*  "igas*  ***^  '7gas-llq^ 
are  also  needed  to  properly  predict  neat  transfer 
rates  and  system  eigenvalues,  that  Is,  the  sinking 
velocity  (vg  or  -v^)  and  the  propellant  surface 
temperature  (Tp,).  The  heat  transfer  coefficients 
tor  qxiq-solld>  ^llq*  ^gas  adopted  In 
the  combined  form  of  forced  convection  and 
conduction.  In  the  early  stage  of  time  period  II 
(described  In  physical  events)  or  before 
propellant  Ignition,  the  radial  velocity  In  Region 
1  of  the  foam  layer  can  be  quite  high.  Therefore, 
the  forced  convection  effect  dominates  heat 


transfer  among  the  spall  particle,  foam  layer,  and 
propellant.  However,  In  the  quenched  cases,  the 
conduction  heat  transfer  Is  not  the  main  mechanism 
during  most  of  time  period  II.  Some  of  the  heat 
transfer  correlations  used  are  given  In  Ref.  12. 

In  consideration  of  various  fwnss  of  heat  transfer 
correlation.  It  Is  noted  that  the  heat  transfer 
measurements  available  for  foam  layers  are  very 
scarce.  In  fact,  heat  transfer  rates  of  foams 
generated  by  liquefaction  and  pyrolysis  of  solid 
propellant  are  completely  unavailable.  This  Is  an 
area  requiring  further  research.  The  heat  fluxes 
mentioned  above  can  be  expressed  as  follows; 


^gas 

*  h 

(T-_ 

-  T,) 

conv 

SB 

g 

•  •1 

'’llq 

*'eonv 

^^SB 

■  V 

(15) 

(16) 


’’llq-solld 


where  T_ 


*'oonv 


'(1*T)p'c 


'pi  * 


The  heat  flux  at  propellant  surface  (qpp) 
Is  calculated  from  the  following  equation,  as 
described  In  Ref.  12. 


(17) 

(18) 


-  3lp  (- 


ps 


-  T, 


pi 


pr 


(19) 


where  the  coefficient  3*p  results  from  a  third- 
order  polynomial  approximation  for  the  temperature 
profile.  The  time  variation  of  the  thermal  wave 
penetration  depth,  ip,.,  is  determined  from 


2M  a 


pr 


8  V  4 


m  pr 


24. 


pi 


)  dt 


(20) 


After  4pp  la  solved  from  Eq.  (20),  the  propellant 
surface  temperature  can  be  determined  from  Eq. 
(11)  through  the  use  of  Eqs.  (17)  and  (19). 


The  sinking  velocity  can  then  be  calculated 
by 

.1,  .« 

V  .  (21) 

*  ^pr^melt 

In  addition  to  the  above  correlations  and 
data  provided  In  Table  1,  a  set  of  input  data  was 
prepared  and  listed  In  Table  2.  The  thermodynamic 
and  transport  properties  were  obtained  from  open 
literature  [IH-IPL 


Table  2.  Input  Data  Used  In  the  HFCI  Simulation 


Region  Parameter  Value  Unit 


Spall 

'•s 

0.3175  X  10"2  or 

0 

0.4851  X  10"2 

m 

Particle 

I'S 

0.3937  X  10"2  or 

a 

0.98  X  10"2 

m 

16.246 

W/m-K 

Cs 

5.02  X  102 

J/kg-K 

Ps 

8.0  X  io3 

kg/ffl3 

2.094  X  10"’ 

H/m-k 

Pi 

1.60  X  103 

kg/m3 

Cl 

1.465  X  103 

J/kg-K 

Foam 

“1 

0.832  X  10"3 

kg/m-s 

Layer 

Ea,tg 

1.597  X  105 

J/mole 

*2 

1.98  X  lo’" 

1/s 

iH* 

-6.695  X  105 

J/kg 

f.t 

AH* 

-1.912  X  10* 

J/kg 

X 

*g 

36.29  X  10"3 

W/m-K 

=pg 

1.465  X  103 

J/kg-K 

n 

0.80 

•••• 

Cpr 

3.852  X  102 
♦  2.598  •  T(K) 

J/kg-k 

Ppr 

1.678  X  103 

kg/ m3 

LOVA 

Ipr 

2.094  X  10*’ 

W/n-K 

Propellant 

"*f.pr 

-7.579  X  105 

J/kg 

^a.Pl 

2.9  X  105 

J/mols 

*1 

1.31  X  1031 

1/s 

Overall  Structure  of  Theoretical  Model 

The  theoretical  model  conaista  of  governing 
equatlona  and  their  aaaoolated  Initial  and 
boundary  eondltlona  for  the  hot  particle,  the 
propellant,  and  the  foam  layer.  Governing 
equatlona  for  the  hot  particle  and  the  propellant 
are  tranatent  heat  conduction  equatlona  written  in 
two-dlmenalonal  cylindrical  form.  Theae  equatlona 
were  recaat  Into  ordinary  differential  equatlona 
Dy  Integral  methoda.  In  theae  two  reglona, 
Inatantaneoua  temperature  dlatrlbutlona  are  aolved 
from  the  governing  equatlona  coupled  to  the  foam 
layer  through  flux  balancea  at  their  boundariea 
{aee  Fig.  5).  The  aolutlon  of  major  unknowna  In 
the  foam  layer  la  alao  delineated  In  the  aame 
figure. 

DISCUSSION  OF  RESULTS 

In  the  following,  one  aet  of  HFCI  almulatlon 
reaulta  la  compared  with  experimental  reaulta,  an.' 
two  aeta  of  parametric  runa  are  preaented  and 


dlacusaed  In  detail.  To  atudy  the  Initial 
temperature  effect  of  the  apall  particle,  two 
different  temperatures  (768  and  1000  K)  were 
considered  for  the  same  spall  particle  size  (Lg  • 
0.3937  cm,  Rg  •  0.3175  cm).  The  calculated 
results  from  the  ease  with  an  Initial  temperature 
of  766  K  are  compared  with  experimental  data 
obtained  by  Miller  and  Cohen  [13].  The  third  set 
of  results  was  obtained  for  studying  the  mass 
effect  on  the  HFCI  process  using  a  large  spall 
particle  (Lg  -  0.98  cm,  rg  •  0.4851  cm}. 

Figure  6  shows  a  set  of  propellant  subsurface 
temperature  profiles  at  different  times  before  the 
spall  particle  begins  to  sink  Into  the  LOVA 
propellant.  As  one  can  see  from  the  temperature 
profile  variation,  the  thermal  wave  penetration 
depth  becomes  deeper  as  time  Increases.  However, 
the  surface  temperature  remains  approximately  at 
the  735  K  range. 

The  axial  temperature  variation  at  the  center 
line  of  the  spall  particle  is  shown  In  Fig.  7  for 
different  times  in  time  period  I.  This  figure 
Indicates  the  decay  of  maximum  temperature  at  a 
distance  far  from  the  spall  particle  base.  The 
temperature  gradient  at  the  base  of  the  spall 
drops  from  a  large  value  upon  Initial  contact,  to 
a  much  lower  level  at  1.3  s.  At  the  end  of  time 
period  I,  the  axial  temperature  distribution 
becomes  quite  uniform. 

In  time  period  II,  the- spall  particle  sinks 
Into  the  LOVA  propellant.  The  calculated  sinking 
distance  versus  time  is  compared  with  measured 
data  In  Fig.  8.  The  calculated  sinking  velocity 
Is  also  shown  in  this  figure.  Note  that  the 
sinking  velocity  (the  slope  of  the  curve) 

Increases  rapidly  at  the  very  beginning,  starts  to 
decrease  at  t  ■  0.75  s,  and  approaches  zero  at  the 
end  of  the  calculation  (t-34  s).  The  final 
calculated  sinking  distance  Is  about  15S  larger 
than  the  measured  one.  This  discrepancy  could  be 
caused  by  errors  In  experimental  measurement  as 
well  as  any  Inappropriate  of  heat  transfer 
correlation  used  in  the  foam  layer  of  the 
theoretical  model.  In  general,  this  comparison 
shows  a  reasonable  agreement  between  calculated 
and  measured  results. 

The  calculated  time  variation  of  the  void 
fraction  In  the  foam  layer  (id)]  la  shown  In  Fig. 
9.  When  the  spall  particle  starts  to  sink  Into 
the  LOVA  propellant,  p  increases  drastically  due 
to  Intensive  heating, -and  then  slows  down 
significantly  as  time  Increases  at  the  end  of 
computation;  only  72*  of  the  space  in  the  foam 
layer  Is  occupied  by  the  gas.  The  calculated  time 
variations  of  center  and  base  temperatures  of  the 
spall  particle  are  compared  wltlt  thermocouple 
measurements  at  corresponding  locations  In  Figs. 

10  and  11,  respectively.  Comparison  of  the 
predicted  base  temperature  Is  better  than 
comparison  of  the  predicted  center  temperature 
with  experimental  data.  Both  of  these 
temperatures  decline  monotonlcally  as  the  particle 
begins  to  sink  into  the  LOVA  propellant  and 
approach  asymptotic  levels  near  the  end  of 
computation. 

The  calculated  time  variation  of  the 
instantaneous  subsurface  temperature  of  the  first 
thermocouple  location  (6  mm  below  the  Initial 
propellant  surface)  compares  well  with  measured 
data,  as  shown  in  Fig.  12.  Like  the  measured 
data,  the  calculated  subsurface  temperature  at  the 
second  thermocouple  location  (12  ram  below  the 
surface)  shows  no  reponse  at  the  end  of 
calculation  (t»3U  s). 


TEMPERATURE  OF  LOVA  PROPELLANT, 
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Instantaneous 
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Distribution 
is  solved  froo 
the  transient 
heat  conduction 
equation. 
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mass  loss 
energy  due  to 
.^ransfer  gasification 


Foam  Layer  X  / _ 

Instantaneous  values  of  porosity,  bulk  temperatures  of  gas 
and  liquid  phases,  radial  velocity  of  foam  layer,  density  of 
the  gas  phase,  and  sinking  velocity  are  solved  from  a  set  of  mass 
and  energy  conservation  equations  and  flux  balances  at  boundaries. 
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Instantaneous  temperature  distribution  la  solved  from  the 
transient  heat  conduction  equation. 


Figure  S.  Block  Diagram  Showing  the  Mathematical  Formulation 

and  Coupling  Relationship  Between  Hot  Particle,  Foam 
Layer,  and  Solid  Propellant 
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Fig.  6  Calculated  Temperature  Distributions  in 
the  LOVA  Propellant  at  Various  Times 


AXIAL  distance  FROM  SPALL  PARTICLE  SASE.I,min 

Fig.  7  Calculated  Axial  Temperature 

Distributions  Along  the  Centerline  of 
the  Spall  Particle  at  Various  Times 


VOID  FRACTIOM  IN  THE  FOAM  LAVER 
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Calculated  Time  Variation  of  the  Void 
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Comparison  of  Calculated  Temperature- 
Time  Trace  at  the  Center  of  tne  Spall 
Particle  with  Experimental  Data 
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Fig.  11  Comparison  of  the  Tenperature-Tlne  Trace 
at  the  Bottom  Surface  of  the  Spall 
Particle  with  Experimental  Data 
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Tiaa  Tracaa  at  th«  First  anc  Second 
Theraocouple  Locations  (6  and  12  aa 
Below  tne  Propellant  Surface)  with 
Experlaental  Data 

Conparlson  of  bulk  temperatures  of  gas  and 
liquid  phases  in  the  foam  layer  Is  shown  In  Fig. 
13.  During  time  period  11,  the  gas-phase 
temperature  is  always  higher  than  that  of  the 
liquid.  The  temperature  difference  decreases 
monotonloally  and  reaches  a  small  difference  near 
the  end  of  computation. 

In  addition  to  the  quenched  case  discussed 
above,  two  separate  HFCl  tests  were  conducted  by 
Miller  and  Cohen  [13],  using  different  spall 
particle  sizes  and  Initial  temperature.  The  small 
cylindrical  particle  has  a  radius  of  0.3175  cm  and 
a  length  of  .3937  cm,  with  an  Initial  temperature 
of  986  K  at  the  time  of  contact  with  the  LOVA 
propellant.  The  large  cylindrical  particle  has  a 
radius  of  O.USSI  cm  and  a  length  of  0.98  cm,  with 
an  Initial  temperature  of  951  K  at  the  time  of 
contact  with  the  LOVA  propellant.  In  both  cases, 
the  propellant  Ignited.  The  pyrolyzed  gases  are 
not  luminous  near  the  base  region  of  the  spall 
particle.  The  observed  luminous  flame  Is  above 
the  Initial  propellant  surface.  The  flame  widens 
abruptly  at  5.2  and  1.2  s,  coresponding  to  the 
small  and  large  particles,  respectively,  widening 
of  the  flame  could  be  considered  as  onset  of 
sustained  Ignition  from  the  experimental  point  of 
view. 

To  validate  the  theoretical  model  with  the 
ignited  cases,  both  tests  were  simulated.  Results 
are  discussed  below.  Figure  lA  shows  the 
comparison  of  calculated  trajectories  of  spall 
particles  with  experimental  data.  The  small 


Pig.  13  Calculated  Temperature-Time  Traces  of 
Gas  and  Liquid  Phases  In  the  Foam  Layer 


Fig.  1A  Corapartson  of  Calculated  Time  Variations 
of  Sinking  Distances  with  Experimental 

Data 
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parclcle  has  a  slightly  higher  Initial  temperature 
than  the  large  one,  and  hence  sinks  Into  the  LOVA 
propellant  at  a  slightly  higher  velocity  (see  Fig. 
15).  At  a  later  time,  the  velocity  of  the  large 
particle  Is  higher  since  It  contains  more  thermal 
energy.  The  calculated  trajectory  shown  In  Fig. 

14  overpredicts  the  experimental  data  by  25S  at 
the  end  of  the  run.  However,  the  agreement  la 
still  acceptable  for  most  of  the  event.  For  these 
two  cases,  the  particles  sink  oohtlnously  as  the 
propellant  reaches  sustained  Ignition  conditions. 
The  calculated  sinking  velocity  variations  with 
respect  to  time,  shown  In  Fig.  15,  exhibit  spikes 
during  the  Initial  Interval  when  the  particles  are 
extremely  hot. 

Figure  16  shows  the  comparison  of  calculated 
temperature-time  traces  at  the  base  of  the  spall 
particle  with  experimental  data.  During  the  early 
phase,  the  base  temperature  for  small  particles  la 
higher  than  that  for  large  particles;  the  trend 
reverses  at  the  later  phase  of  Time  Period  II 
(when  spall  particles  sink  Into  the  LOVA 
propellant).  Calculated  results  are  In  good 
agreement  with  experimental  data  In  terns  of  the 
trend  mentioned  above,  as  well  as  the  magnitude 
and  slope  of  these  traces.  It  Is  Interesting  to 
note  that  the  average  temperature  at  the  spall 
particle  base  drops  continuously  even  after  or.set 
of  sustained  ignition;  this  la  believed  to  be  the 
result  of  Incomplete  combustion  In  the  foam  layer. 
Caseous  pyrolysis  prouuots  In  the  bubbles  take 
time  to  reach  the  fully  reacted  state.  Heat 
release  In  the  foam  tone  la  sufficiently  low  to 


Fig.  15  Calculated  Time  Variations  of  Sinking 
Velocity 


Fig.  16  Comparison  of  Calculated  Temperature- 
Time  Traces  at  the  Base  of  the  Spall 
Particle  with  Experimental  Data 


produce  a  relatively  cool  layer  of  foam  material 
which  reduces  the  energy  of  the  spall  particle. 

The  above  effect  can  be  seen  from  the  plot  of  the 
calculated  bulk  temperatures  of  the  gas  and  liquid 
phases  In  the  foam  layer  (see  Fig.  17).  The 
temperatures  of  both  phases  are  lower  than  the 
average  temperature  of  the  spall  particle  base. 

The  calculated  void  fraction  variations  with 
respect  to  time  Is  shown  In  Fig.  18.  After  the 
Initial  Intensive  heating  In  time  period  II,  the 
void  fraction  reduces  slightly  with  respect  to 
time. 

Comparisons  of  calculated  temperature-time 
traces  at  the  center  of  spall  particles  with 
experimental  data  Is  shown  in  Fig.  19.  Agreement 
Is  quite  reasonable,  and  the  trend  Is  Identical  to 
that  described  for  spall  particle  base.  Before 
completion  of  the  computation  of  10  s  for  the  time 
period  II.  the  center  temperatures  are  always 
higher  than  those  at  the  spall  particle  base. 

Based  upon  these  comparisons.  It  Is  quite 
obvious  that  the  HFCI  model  Is  able  to  predict  the 
quenched  and  Ignited  cases  using  different  spall 
particle  sizes.  Also,  by  comparing  results  of  the 
small  particle  (L  ■  0.3937  cm,  r^  •  0.3175  cm)  at 
an  Initial  temperature  of  986  K  with  those  for  the 
same  sized  particle  at  1000  K,  one  can  observe  the 
strong  difference  In  solution  from  a  quenched  case 
to  a  runaway  Ignition.  According  to  the  test  data 
of  Hiller  and  Cohen  [13],  the  small  particle  can 
Ignite  the  LOVA  propellant  at  986  X,  which  is  In 
agreement  with  the  prediction  of  ignition  at 
986  K. 


TEMPERATURE  OF  THE  LIQUID  PHASE  BULK  TEMPERATURE  OF  THE  GAS  PHASE 

IN  THE  FOAM  LAYER.  T,.  K  .  'N  THE  FOAM  LAYER,  T,.  K 
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SUMMARY  AMD  CONCLUSIONS 

1)  The  two-step  chemical  kinetic  mechanism, 
proposed  by  Hiller  et  al.  [8],  based  upon  DSC 
measurements,  was  Incorporated  Into  the 
theoretical  model  [12]  for  simulating  hot 
fragsent  conductive  Ignition  processes  of 
LOVA  propellants. 

2)  Calculated  results  compared  well  with 
experimental  data  In  temperature-time  traces 
of  spall  particle  and  propellant. 

3)  Reasonable  agreement  was  achieved  between 
theoretical  prediction  and  experimental 
measurements  In  terms  of  the  terminal  position 
and  particle  trajectory  of  the  quenched  spall 
particle  which  Is  partially  submerged  in  the 
LOVA  propellant.  The  predicted  trajectory  of 
the  spall  particle  also  agrees  quite  well  with 
measured  data  for  Ignited  cases. 

V)  The  controlling  mechanism  In  the  HFCI  process 
Is  the  competition  between  the  rate  of  heat 
release  from  exothermic  reactions  Introduced 
by  the  spall  particle  and  the  rate  of  heat 
loss  to  the  surrounding. 

5)  The  present  two-step  reaction  model  gave  a 
reasonable  numerical  solution  In  HFCI 
simulation.  However,  a  more  complicated 
reaction  model  may  be  needed  for  various  types 
of  LOVA  propellants.  Also,  additional  heat 
transfer  correlations  in  the  foam  zone  of  the 
pyrolysis  product  of  various  propellants 
should  be  developed  for  more  accurate 
predictions  of  HFCI  processes. 

7)  From  HFCI  point  of  view.  It  Is  highly 

desirable  to  have  LOVA  propellants  with  cool 
pyrolysis  products  In  the  foam  layer  which  can 
effectively  quench  the  not  spall  particle. 
Also,  It  la  desirable  to  have  LOVA  propellants 
which  require  extended  distances  for 
development  of  flames. 
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